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i , The discovery of the top quark in 1995, by the CDF and D0 collabora- 

(D [ tions at the Fermilab Tevatron, marked the dawn of a new era in particle 

' physics. Since then, enormous efforts have been made to study the properties 

of this remarkable particle, especially its mass and production cross section. 
In this article, we review the status of top quark physics as studied by the 
two collaborations using the pp collider data at i/s = 1.8 TeV. The combined 
measurement of the top quark mass, mt = 173.8±5.0 GeV/c 2 , makes it known 
to a fractional precision better than any other quark mass. The production 
cross sections are measured as a t j = 7.61^'g pb by CDF and a t j = 5.5± 1.8 pb 
by D0. Further investigations of it decays and future prospects are briefly 
discussed. 
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1. INTRODUCTION 

The discovery0~i of the top quark in 1995 was a major triumph of the Standard Model 
of particle physics.S It was the culmination of nearly two decades of intense research at 
accelerators around the world. The direct measurement of a large mass for the top quark, 
by far the heaviest fundamental particle known, has caused much excitement. That the mass 
is close to the electroweak scale suggests the tantalizing possibility that the top quark may 
play a role in the breaking of electroweak symmetry and therefore in the origin of fermion 
masses. The top quark mass is one of the most important parameters of the Standard Model. 

Since the discovery, the CDF and D0 collaborations have collected more data and per- 
formed detailed studies. They have refined particle identification techniques and adopted 
innovative analysis methods, resulting in precise measurements of the top quark mass and 
the it production cross section. There exist several excellent reviews! of the work that led to 
the discovery and the work done shortly thereafter. Although we shall touch upon some of 
the highlights of that exciting time, the focus of this review is the current status of top quark 
physics resulting from the recent measurements from the Tevatron experiments. 

In the remainder of this section, we give a sketch of the Standard Model, followed by 
discussions of the arguments and evidence for the existence of the top quark that predate 
the discovery, the indirect measurements of the top quark mass from electroweak data, and 
the significance of the heavy top quark in taking us beyond the Standard Model. In Sec. [2|, 
we outline the top quark production mechanisms in pp collisions, the decay signatures of 
it pairs and the Monte Carlo modeling of it events. The ingredients involved in making 
and detecting the top quarks, such as the Tevatron collider complex, the detectors, particle 
identification techniques, and the characteristics of the signal and background are discussed 
in Sec. |[ The saga of the early searches for the top quark and of its discovery at the Tevatron 
is summarized in Sec. f|. Sections || and || describe, respectively, the measurements of the it 
production cross section and the top quark mass by the CDF and D0 collaborations. Several 
other studies that have been made using the present it event samples are summarized in 
Sec. [F[ Finally, in Sec. ||, we discuss the prospects for the next collider run, scheduled to 
begin in the year 2000. 



1.1. Synopsis of the Standard Model 

The Standard Model (SM), the prevailing theory of matter and forces, has been in place 
for over two decades. The particles of matter are spin-1/2 quarks (q) and leptons (£), which 
seem to be elementary, at least down to 10~ 18 meters. There are six "flavors" of quarks, and 
likewise of leptons, grouped in pairs into three generations. They interact via the exchange 
of spin-1 gauge bosons: eight massless gluons, the massless photon, and the massive 
and Z° bosons. The top quark was the important missing piece in the fermion sector of the 
Standard Model. The building blocks of the Standard Model are shown in Table [1]. 

A vital part of the Standard Model that awaits experimental evidence is the "Higgs 
mechanism." The Standard Model is based on the gauge group SU(3)c x SU(2)l x U(1)y 
and accommodates electroweak and flavor symmetry breaking by introducing a weak-isospin 
doublet of fundamental scalar fields $ = y, J with the potential function 
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Table 1. Particles of the Standard ModelJl 





symbol 


name 


mass (MeV/c 2 ) 


charge (e) 


On ark d 

(spin = 1/2) 


u 
d 


11Tt 

down 


ft 

« 10 


2/3 

-1/3 
/ u 




Q 


rbarm 

v. .X ±( (1 X. 1 J. 


« 1500 


2/3 




s 


strantre 


« 200 


-1/3 




t 

b 


top 
bottom 


173 8 GeV/c 2 
w 4500 


2/3 
-1/3 


Leptons 


^ e 


electron neutrino 


< 10 eV 


o 


(spin=l/2) 


e 


electron 


0.511 


-1 




V,, 
A* 


muon neutrino 


< 0.17 


o 




U 


muon 


105.7 


-1 




V T 


tau neutrino 


< 24 







T 


tau 


1777 


-1 


Gauge bosons 


7 


photon 








(spin = 1) 


W 

z 


W 
Z 


80.3 GeV/c 2 
91.2 GeV/c 2 


1 






9 


gluon 








Higgs 


H 


Higgs 


? 


? 



V($ f $) = /i 2 ($ t $) + iak^ 1- *) 2 , (1) 

where A is the self coupling of the scalar field. With /z 2 chosen to be negative, the electroweak 
symmetry is spontaneously broken (that is, the vacuum state fails to display the symmetry 
of the theory) when the scalar field is expanded about its (non-zero) vacuum expectation 
value v = \J / A = (Gp\^2)~^ = 246 GeV (referred to as the electroweak scale). The 
spontaneous breaking of electroweak symmetry endows the W ± and Z° bosons with masses, 
= na/CpV^sm 2 9 W and Mf = M^j cos 2 9w and also gives rise to a spin-0 (scalar) 
particle called the Higgs boson. (Here, a is the fine structure constant, Gp is the Fermi 
(weak) coupling constant, and 6w is the weak angle.) Each quark and lepton / has its own 
Yukawa coupling to the Higgs boson Gj and thus acquires a mass rrif — Gfv/ y/2. 

The Standard Model has been extremely successful, so farP It has withstood scores of 
very stringent experimental tests in a variety of high-energy interactions. No significant 
discrepancies between experimental data and the Standard Model have yet been found. 
But several critical issues remain unresolved. The inclusion of the Higgs mechanism is 
artificial: there is no explanation for the form of the Higgs potential, and therefore for 
neither electroweak symmetry breaking nor the breaking of flavor symmetry. Hence, physics 
beyond the Standard Model seems inevitable,i and it is entirely plausible that the top quark 
might be our window to that new physics. 

1.2. Why Must the Top Quark Exist? 

Long before the top quark was observed, there were compelling arguments for its ex- 
istence. The renormalizability of the Standard Model requires the cancellation of triangle 
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z 

Fig. 1. An example of a triangle diagram that gives rise to a chiral anomaly in the Standard 
Model. 



anomalies — a problem that arises from the interaction of three gauge bosons via a closed 
loop of fermions as shown in Fig. |I[ It turns out that the fermion contributions within each 
generation cancel if the electric charges of all left-handed fermions sum to zero: 



■1 + 3 x 



0. 



(2) 



The factor 3 is the number of color charges for each quark flavor. For this to work for the 
third generation, the top quark with Q = 2/3 must exist. 

There is ample indirect experimental evidence for the existence of the top quark. The 
experimental limits on flavor changing neutral current (FCNC) decays of the 6-quarki such 
as b — > s£_L£~ and the absence of large tree level (lowest order) B®B® mixing at the T(45) 
resonance^O rule out the hypothesis of an isosinglet 6-quark. In other words, the 6-quark 
must be a member of a left-handed weak isospin doublet. 

The most compelling experimental evidence comes from the wealth of data accumulated 
at e + e~ colliders in recent years, particularly the detailed studies of the Zbb vertex near the 
Z resonance. These studies have yielded a measurement of the isospin of the 6-quark. The 
Z boson is coupled to the 6-quarks (as^ell as to other quarks) through vector and axial 
vector charges (v b and a b ) with strengtl 




G F M\ 
\ 2V2 



(3) 



where v b and a b are given by 



v b = 2[l£(b)+li 



— Ae b sin 2 8 W , 



and 



(4) 



a b = 2[l£(b) + I«(b)}. (5) 

Here, I%(b) and I3 (b) are the third components of the isospin for the left-handed and 
right-handed 6-quark fields. The electric charge of the 6-quark, e b = —1/3, has been well 
established from the T leptonic width as measured by the DORIS e + e~ experiments.^ The 
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Born approximation in the limit of a zero mass 6-quark gives for the partial Z boson decay 
rate 

T bl ^T{Z^bb)^^-^{vl + al). (6) 

The partial width T bb is expected to be thirteen times smaller if l^(b)=0.0. The LEP 
measurement of the ratio of this partial width to the full hadronic decay width, R b = 
Tfe/Thad = 0.2170 ± 0.0009, is in excellent agreement with the SM expectations (including 
the effects of the top quark) of 0.2158,E-3 ruling out I^(b)=Q.Q. In addition, the forward- 
backward asymmetry in e + e~ — > bb at the Z resonance, 

_ 3a e v e a b v b 

~ (v* + al)(vl + aiy (?) 

is sensitive to the relative size of the vector and axial vector couplings of the Zbb vertex. 
The sign ambiguity for the two contributions can be resolved by the Aps measurements 
from low energy experiments that are sensitive to the interference between neutral current 
and electromagnetic amplitudes. So, from the measurements of T bb and Apb at LEP, SLC, 
and the low energy experiments (PEP, PETRA and TRISTAN), one obtains^ 

Ii(b) = -0.490j™& (8) 
I*(b) = -0.028 ± 0.056, 

for the third component of the isospin of the 6-quark. This implies that the 6-quark must 
have a weak isospin partner, i.e., the top quark, with I$(t) = +1/2. 



1.3. Indirect Constraints on the Top Quark Mass 

An upper bound on the top quark mass can be obtained by requiring that partial wave 
unitarity be respected at tree level in the reactions it — > W + W~ , ZZ, HZ, and HH. This 
leads to a condition on the top quark mass, which sets the scale m t {Gp\/2)~^ of Hit Yukawa 
couplings and a constraint m t < 500 GeV/c 2 .0 

Since virtual top quarks are involved in higher order electroweak processes, tighter con- 
straints on the top quark mass can be obtained from precision electroweak measurements. 
The higher-order (radiative) corrections to many electroweak variables depend on the masses 
of the top quark and Higgs boson via loop diagrams such as those shown in Fig. |. 

At one loop, for example, the p parameter, 



P 



M 2 W 



MlU-sm 2 6i 



1 + Ar, 



l ZK L ~ 0111 °W) 

which relates the W and Z boson masses and the weak angle, gets a radiative correction 



(9) 



Ar 



3Gp 2 \f^2Gp 
- r m t + 



IQty 



2 



11 (M 2 H 



+ 



+ 



(10) 
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H 




Fig. 2. Self-coupling loops contributing higher order quantum corrections. 

which is quadratic in the top quark mass. Note, however, that the dependence on the mass 
of the Higgs boson Mh is only logarithmic. Therefore, the top quark mass, especially if 
large, is the dominant parameter in corrections to electroweak processes. This relation was 
used to set early constraints on the mass of the top quark;0 for example, ignoring the Higgs 
contribution, Ar < 0.02(0.01) implies that m t < 250(180) GeV/c 2 . Additional constraints 
can be derived from the large body of precision electroweak data. Taking as inputs the 
extremely high precision measurement of the Z boson mass, sin 2 9w from the Z boson's 
decay rate, the forward-backward asymmetry and left-right polarization measurements in 
Z boson decay, and the sin 2 9w measurement from the vN scattering, a to the Stan- 
dard Model predictions with m t and Mh as free parameters yields m t = 157ig° GeV/c 2 
and Mh = GeV/c 2 . When the Mw measurements from the Tevatron and LEP are 

included, the resulting mass@ of the top quark is m t = 181. 31^2 jZif;| GeV/c 2 , where the 
first uncertainty is, as before, the statistical error from the fit and the second uncertainty 
reflects the assumed variation of Mh in the range 70-1000 GeV/c 2 . These precision data, 
along with the direct measurements of the masses of the top quark and the W boson, provide 
an indirect measurement of the mass of the Higgs boson. 



1.4. Significance of a Heavy Top Quark 

The top quark is the heaviest elementary particle yet discovered. Its mass, of the same or- 
der as the electroweak scale (~ v/y/2), is about twice that of the W and Z bosons and about 
40 times larger than its isospin partner, the 6-quark. Because of its large Yukawa coupling to 
the Higgs boson (G t ~ 1), and hence to the mechanism of electroweak symmetry breaking, 
the top quark may have unique dynamics. Its mass has already set severe constraints on 
extensions to the Standard Model, including any new theories of strong interactions, leading 
to the development of top condensation, topcolor, and related ideas. 

But the most intriguing observation of all is that in supersymmetric models with grand 
unification, a large top quark mass will automatically break electroweak symmetry in the 
required manner .1121 At the grand unification scale Mqut, well above the weak scale Mw, all 
the supersymmetric scalars (the squarks and sleptons, denoted /) will have the same mass: 



9 



M?(M GUT ) = ml (11) 

The two Higgs scalars will also have the same mass: 

Ml H {M GVT ) =m 2 + fi 2 . (12) 

As one moves to smaller energy scales, these masses evolve according to the renormalization 
group equations. For the mass of the Higgs scalar at a scale Q, one finds: 

M 2 (M GUT ) - Ml(Q) oc m\ In ■ (13) 

For a sufficiently large top quark mass {m t ~ 175 GeV/c 2 ), it is therefore possible that 
M%(Mw) < at the weak scale, which is required to break electroweak symmetry. The 
squark masses evolve in a similar manner, but with a smaller proportionality constant. 
Thus, one can avoid breaking the SU(3) color symmetry. 



2. TOP QUARK PRODUCTION AND DECAYS 

The dominant production mechanism for top quarks at a hadron collider is pair produc- 
tion (that is, tt). They can also be produced singly, but with a rate calculated to be about 
half that for pair production. The main characteristics of these processes are as follows: 

1. Pair production of top quarks through the quantum chromodynamic (QCD) pro- 
cesses qq — > it and gg — > tt. (See Fig. |3|.) At the Tevatron, the relative contributions of 
these two processes are about 90% and 10%, respectively. There are also contributions 
with intermediate photons or Z bosons, but they are much smaller and can safely be 
ignored. 

2. Drell-Yan production of single top quarks through pp — > W* + X — > tb + X. 
(See Fig. ^.) This process would have been dominant at yfs = 630 GeV if m t were 
less than (Mw — mb) , which is why earlier top quark searches at CERN0 were based 
on detecting this mode (see Sec. |j). However, it still contributes significantly to the 
inclusive top quark production cross section for m t ~ 175 GeV/c 2 . 

3. Single top quark production via ty-gluon fusion. (See Fig. |5|.) Photon-gluon 
and Z-gluon processes are also allowed, but, again, the rates are very small. 

In the rest of this section we will discuss QCD pair production of top quarks. The single 
top quark production modes will be discussed in Sec. 



2.1. QCD Pair Production of Top Quarks 

Several reviewsilil cover recent developments in the calculation of heavy-quark produc- 
tion cross sections. Here, we will summarize the main issues involved in the calculation of 
the tt production cross section in perturbative QCD. 
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The total production cross section for pp — > tt, where the proton and antiproton each 
have momentum P (= y/s/2), can be factorized in the standard way@ 

a{jjp ~+tt) = ^2j dx a dx b fP(x a , fi 2 )ff(x b , /j, 2 )a(ab -> tt; s, /j 2 , m t ), (14) 

a,b 

where the summation indices a and b run over light quarks and gluons. This formula 
expresses the total cross section in terms of the parton-parton processes ab —>■ it, where 
a and b are partons contained in the initial proton and antiproton carrying momentum 
fractions of x a and x b , respectively. The parton distribution functions and fjf are the 
probability densities of finding a parton with a given momentum fraction in a proton or 
antiproton, and a is the subprocess cross section at a parton-parton center-of-mass energy 
of s = 4x a x b P 2 = x a x b s. The renormalization and factorization scales, here chosen to be 
the same value fi, are arbitrary parameters. The first is introduced by the renormalization 
procedure, and the second by the splitting of the total cross section into perturbative (a) 
and nonperturbative (f p , p) parts. The dependence of observables on fi is an artifact of 
truncating the perturbation expansion at finite order; if the calculations could be carried 
out to all orders, the dependence on fi would vanish. For tt production, one usually takes 
\i ~ m t . Theorists typically estimate the error introduced by truncating the series by varying 
\i within some (arbitrary) range, such as m t /2 < fi < 2m t . However, it should be recognized 
that neither the choice of /i, nor the range over which it is allowed to vary, have physical 
significance. 

The first calculations of the parton-parton cross section a were to 0(a 2 ), that is, to 
leading order (LO)." At the Tevatron, the qq — > tt process dominates, contributing 90% of 
the cross section, while the gg process contributes the other 10%. (The difference between 
the strengths of these two subprocesses arises mainly from differences in the parton distri- 
bution functions for quarks and gluons, rather than from differences in the parton-parton 
cross sections.) Subsequently, several groups calculated the complete 0(af) next-to- leading 
order (NLO) cross sectionES The NLO cross section is about 30% higher than the LO cross 
section; the estimated uncertainty (from the sensitivity to variations in \i and changes in 
the parton distribution functions) is typically 10-20%. 

In a regime where perturbation theory is valid, the NLO contribution should be small 
compared to the LO terms. However, for top quark production at the Tevatron, the NLO 
contribution is worryingly large: for the gg process it is about 70% of the size of the LO 
terms.il (The situation is better for the qq process, where the NLO contribution is about 
20% that of LO.) The large difference between the LO and NLO calculations is mainly due 
to processes involving the emission of soft initial state gluons. Fortunately, it is possible, 
through a technique called resummation, to calculate the sums of the dominant logarithms 
from soft gluon emission to all orders inperturbation theory. This was first carried out by 
Laenen, Smith, and van Neerven (LSvN).o A difficulty arises, however, in that the resummed 
gluon series is divergent due to nonperturbative effects as a s becomes large. LSvN solved 
this problem by introducing a new scale /x 3> Aq CD which is used as a cutoff to remove this 
divergence. They predict an increase in the cross section by 10% over the NLO prediction; 
the uncertainty, however, is relatively large (> 10%), due to the dependence on /i . 

More recently, two other groups have performed this calculation using methods which 
avoid the need for an arbitrary cutoff. Berger and Contopanagos (BC),&~§1 use the technique 
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Table 2. Results of several different pp — ► it calculations, for m t = 175 GeV/c 2 and 
-v/s = 1.8 TeV. Note that the LSvN result uses an older set of structure functions, which makes it 
systematically low compared to the other results.il 



Calculation Type Structure Function a t j 



(1) Exact NLollI 


NLO only 


MRSR2tl 
MRSD'il 


4.871^ pb 


(2) LSvN@ 


Resummed 


4.94lg;S Pb 


(3) BC§ 


Resummed 


CTEQ30 


5.52lg;g pb 


(4) BCMNlIl 


Resummed 


MRSR2H 


5-06±8;|| P b 



of principal value resummation£2l They also find an increase of about 10% over the NLO 
prediction, with an estimated uncertainty of about 5%. Bonciani, Catani, Mangano, Nason, 
and Trent adue (BCMNT)Hl use a slightly different scheme to avoid the divergence and 
treat the subleading log terms differently. They find a much smaller enhancement of the 
NLO cross section, on the order of 1%, with estimated uncertainties also of about 5%. A full 
discussion of the differences between these calculations is beyond the scope of this review, 
but the subject has been discussed in detail in the literature.ilSS Note, however, that 
if one is comparing these calculations to the present experimental results, the discrepancy 
between them is of no practical importance, as the difference between them is substantially 
smaller than the uncertainties on the experimental measurements (~ 30%). The results of 
various calculations for m t = 175 GeV/c 2 are compared in Table |2|. A plot of the various 
cross sections is given in Fig. ||. 

It should be appreciated that these cross sections are extremely small - - about ten 
orders of magnitude smaller than the total inelastic pp cross section. Of the five trillion or 
so collisions which have occurred at each of the CDF and D0 interaction regions during 
Run 1, one expects about 600 it pairs to have been produced. It should also be noted that 
all these calculations assume a minimal form of the Standard Model. Certain extensions to 
the SM, such as models with two Higgs doublets or supersymmetry, predict it production 
cross sections which are different from that of the SM by a few percent.il 



2.2. Top Quark Hadronization and Decay 

Within the Standard Model, the dominant decay of a top quark is via t — * Wb, with a 
branching ratio of nearly 100%. The decays t — > Ws and t — > Wd are also allowed, but are 
suppressed by factors of 10~ 3 -10~ 4 by the Cabibbo-Kobayashi-Maskawa (CKM) mixings.^ 
Other decays, such as flavor-changing neutral current (FCNC) decays, are predicted to be 
many orders of magnitude smaller. We will discuss these further in Sec. [7[ 

Owing to the large mass of the top quark, its lifetime is extremely short (r m 1.5 GeV, 
corresponding to r ~ 5 x 10~ 25 s); so short, in fact, that at Tevatron energies it decays 
before it has a chance to hadronize.Ea This implies that a decaying top quark can be treated 
as a free particle. Note, however, that it isn't just the large mass of the top quark which 
gives it its short lifetime, but also the fact that it has a CKM-allowed decay into a 6-quark. 
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Fig. 6. Results of top quark production cross section calculations for ^/s = 1800 GeV. The 
upper group of curves are for it production,§!lH!i!@ and the lower curve is for single top quark 
production£3 Note that the LSvN result uses an older set of structure functions, and is thus 
systematically low compared to the other results.il 
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Table 3. Decay modes for a tt pair and branching fractions. 





(1/9) 


(1/9) 


W -> rz/ T 
(1/9) 


IV — > qq 
(2/3) 


w - 


-> qq 


(2/3) 


12/81 


12/81 


12/81 
It -4- ipts i 


36/81 
fall iets) 


w - 




(1/9) 


2/81 
(er) 


2/81 
(fir) 


1/81 
(tt) 




w - 


"> 


(1/9) 


2/81 
(e/i) 


1/81 
(/i/i) 






w - 


-> pz/ e 


(1/9) 


1/81 

(pp) 









A fourth-generation down-type quark b' of mass comparable to the top quark might still 
have a long lifetime if all its decay modes were suppressed by the CKM mixings. 

A it final state contains two t — > Wb decays. The two 6-quarks will form jets, while 
each W boson will decay into either a lepton-neutrino or a quark- ant iquark pair. To a 
good approximation, each possible decay of the W boson is equally probable; however, one 
must remember to count each quark flavor three times, since quarks come in three colors. 
Therefore, the probability for a W boson to decay into each of the three lepton flavors is 
about 1/9, while the probability for it to decay into the two available quark final states is 
about 2/3. 

Since there are two top quarks in each event, and since the W bosons decay independently 
of each other, the events can be classified according to how the W bosons decay (see Table 0). 

• Events in which both W bosons decay leptonically are called dilepton events. Since 
tau leptons are difficult to identify, the particular dilepton channels which have been 
most studied are the pp, /i/i, and p/t channels. These final states have the signature of 
two high-pr leptons, a large imbalance in the total transverse momentum ("missing- 
Et" or $ T ), and two 6-jets. These events are expected to have small backgrounds 
(especially the e/i channel). However, as can be seen from Table |3], they also have small 
branching fractions, with all three of these channels comprising only about 4/81 ~ 
4.9% of tt decays. (There is also a small contribution from the r lepton channels. For 
example, it — ► bbevTv — > bbev evv can contribute to the pp channel.) Dilepton events 
also have the drawback of containing two unobserved neutrinos in the final state, which 
prevents complete reconstruction of the event kinematics. 

Several recent analyses have also considered the pt and fir channels. These will be 
discussed in Sec. [572]. 

• Events in which one W boson decays leptonically and the other decays into quarks are 
called lepton+jets events. Those which have been studied are the p + jets and /t + jets 
channels. They are characterized by a final state containing one high-pr lepton, large 
$ T , and four jets, two of which are 6-jets. Compared to the dilepton channels, the 
lepton+jets channels have a much larger cross section — the branching ratio for each 
is about 4/27 ~ 15%. (Again, there is a small additional contribution from r lepton 
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channels.) The final state contains only one neutrino, so there is sufficient information 
to completely reconstruct the event (once a particular set of assignments of jets to 
the final state partons is assumed). The disadvantage of these channels, however, is 
a large background from inclusive W boson production with associated jets, plus a 
smaller background from QCD jet production. 

• Events in which both W bosons decay into quarks are called all-jets events. This final 
state consists of six jets, of which two are fe-jets, no high-p^ leptons, and small $ T . 
This channel boasts the largest branching ratio (~ 44% of the total). Unfortunately, 
that is more than countered by a huge background from QCD multijet processes. 
Nevertheless, it is still possible to isolate a it signal in this channel. The techniques 
for doing so will be discussed in Sec. [5|. 



2.3. Modeling Top Quark Events 

Accurate modeling of the kinematics of it production and decay is essential for extracting 
reliable information from the data. The most widely used general-purpose model is that 
provided by the HERWIG Monte Carlo program.l§ HERWIG models it production starting 
with the leading-order hard process, choosing the parton momenta according to the weight 
given by the matrix element of the process. Gluon emission from both the initial and 
final states is modeled using leading-log QCD evolution,^ keeping track of the correlations 
induced by color strings between partons. Each top quark is then decayed to a W boson 
and a 6-quark, and partons remaining in the final state are hadronized into jets. Products 
of interactions among the beam remnants, called the underlying event, are also included in 
the model. Detector effects are then added using a model of the detector response to the 
physical objects. 

There have been several studies comparing the predictions of parton shower Monte Carlo 
programs such as HERWIG to the more explicit calculations. Frixione et al.0 compare HER- 
WIG to the full NLO calculation. They distinguish two types of quantities: those which are 
delta functions at leading order, such as Pt{H) and A(p(ti), and those which are nontrivial at 
leading order, such as Pr(t), and m t j. For the latter, nontrivial, set of quantities they 
find good agreement between the Monte Carlo predictions and the full NLO calculation. For 
the other quantities, however, they find significant disagreement in the low-pr region (see 
Fig. ^). They interpret this as a deficiency in the NLO calculation due to the lack of resum- 
mation effects, and conclude that HERWIG is more reliable than the NLO calculation in that 
region. Orr et aI.El compare an 0(a^) calculation of gluon emission in it events to HERWIG; 
they find that HERWIG seems to generate too much final state radiation. On the other hand, 
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Dare a resummed calculation for tt production to a different Monte Carlo 



program, pythiaJ^ and find that the latter generates too little radiation. These discrep- 
ancies are probably not large enough to be important for the current experimental results. 
But they will have to be understood better for the next round of experiments — radiation 
effects are already one of the dominant uncertainties for the present measurements of the 
mass of the top quark, so the precision of the measurement may not improve much until 
these effects are better understood. 
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Fig. 7. Transverse momentum of the tt pair: NLO calculation compared to HERWIG. The inset 
magnifies the low pt portion of the plot on a linear scale. The HERWIG prediction has been scaled 



by a constant factor of K = 1.34. From Ref. 22. 
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3. DETECTING THE TOP QUARK 

In this section, we discuss some of the experimental aspects of detecting the top quark. 
We first summarize the apparatus used: the Tevatron collider and the CDF and D0 de- 
tectors. We then discuss the procedures used by the experiments to identify final state 
objects, such as electrons, muons, and jets. We conclude this section with a summary of the 
distinctive characteristics of the top quark signal and its principal backgrounds. 



3.1. The Accelerator 

The Fermilab pp collider^ in Batavia, Illinois is the world's highest energy particle 
accelerator, with a center-of-mass energy of 1800 GeV. It is the only facility at present 
capable of producing top quarks for direct study. 

A schematic of the accelerator complex is shown in Fig. |j. Negative hydrogen ions are 
first accelerated to 750 keV by an electrostatic Cockroft-Walton accelerator and then further 
boosted to 400 MeV by a 150 m long linear accelerator. At the end of this accelerator, the 
electrons are stripped from the ions and the resulting protons enter the booster. This 75 m 
radius synchrotron accelerates the protons to 8 GeV. From there they are injected into the 
Main Ring. 

The Main Ring is a large (1000 m radius) synchrotron, composed of conventional elec- 
tromagnets, which accelerates protons and antiprotons up to 150 GeV for injection into the 
Tevatron. It also serves as a source of 120 GeV protons for producing antiprotons.0 The an- 
tiprotons are collected from the production target using a lithium lens, momentum-selected 
around 8 GeV, and then directed first into the Debuncher and then into the Accumulator. 
These are two concentric storage rings with radii of about 80 m. There, the antiprotons 
are stochastically cooled^ to reduce their momentum spread. When enough antiprotons 
have been accumulated (stacked), they are extracted into the Main Ring, accelerated, and 
injected into the Tevatron. 

The Tevatron, a synchrotron made from superconducting magnets, is situated just below 
the Main Ring. In collider mode, the Tevatron is filled with six bunches of protons and six 
bunches of antiprotons, circulating in opposite directions. The beams are accelerated to 
the maximum energy of 900 GeV each and brought into collision at the CDF and D0 
experimental areas. The beams are typically kept colliding for about 20 hours, after which 
the machine is emptied and refilled with new batches of protons and antiprotons. The length 
of each bunch is about 50 cm, dictated by the accelerator RF system, giving a luminous 
region at each interaction point which is roughly Gaussian with a longitudinal width of about 
30 cm. This relatively long bunch length produces a degradation of the transverse energies 
at the trigger level because the information about the position of the interaction point is 
not available at the early stages of triggering. At CDF, it also results in a substantial loss 



of acceptance of the vertex detector for 6-quarks (see Sec. |3.3.6.2j ). 

The Main Ring lies mostly in a plane, except at the CDF and D0 experimental areas 
where it is bent into overpasses to allow room for the detectors. The separation between 
the Main Ring and the Tevatron is 5.8 m at CDF and 2.3 m at D0. At CDF, the overpass 
clears the detector. At D0, however, the Main Ring overpass goes through the outer (coarse 
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Table 4. Parameters of the Fermilab Tevatron collider for Run 1. 

Accelerator radius 
Maximum beam energy 
Injection energy 
Peak luminosity 
Number of bunches 
Intensity per bunch 
Crossing angle 
Bunch length 
Transverse beam radius 
Fractional energy spread 
RF frequency 
p stacking rate 
Beam crossing frequency 
Period between crossings 



hadronic) part of the calorimeter. This is unfortunate because during normal collider data- 
taking, the Main Ring is used for antiproton production and losses from the Main Ring 
may deposit energy in the detectors, thereby increasing background. D0 rejects much of 
this background at the trigger level by rejecting triggers that occur during injection into the 
Main Ring, when losses are large. Some triggers are also disabled whenever a Main Ring 
bunch passes through the detector or when losses are registered in scintillation counters 
around the Main Ring. This results in a loss of about 15% of the available livetime. The 
problem of Main Ring contamination is far less severe at CDF; even so, CDF observes 
occasional events with extra energy from the Main Ring. These events are rejected offline. 

The collider was commissioned with a short run in 1985, followed by the first high 
luminosity run in 1988-1989. Only the CDF detector took data during that run, which 
had a peak luminosity of about 2 x 10 30 cm' 2 s _1 (a factor of two greater than the design 
luminosity). The second series of runs took place over the period 1992-1996, during which a 
peak luminosity of 2 x 10 31 cm~ 2 s _1 was achieved. This period was divided into three runs, 
designated Run la, Run lb, and Run lc; the delivered integrated luminosities for these runs 
were about 23 pb _1 , 122 pb _1 , and 17 pb -1 , respectively. The results covered in this review 
are from this running period. Table |] reviews the major parameters of the collider. 

3.2. The CDF and D0 Detectors 

The two collider experiments at the Fermilab Tevatron, CDF§@ and D0,i are illus- 
trated in Figs. [| and [IT]. Both were designed to study high-p T interactions and feature large 
angular coverage and good identification and measurement of electrons, muons, and jets. 
The layouts of both detectors are broadly similar. Moving outwards from the interaction 
point, one first encounters tracking detectors, which measure the trajectories of charged 
particles, then calorimeters, which measure the energies of jets and of electromagnetic show- 
ers, and finally an outer set of tracking chambers, which identify and measure muons that 
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CDF Detector 




INTERACTION REGION 



Fig. 9. The CDF detector, in isometric view (left) and cross sectional view of one quadrant 
(right). The detector is forward-backward symmetric about the interaction region, which is at the 
lower-right corner of the cross sectional view. SVX, VTX, CTC, and CDT are tracking detectors, 
CEM, CHA, WHA, PEM, PHA, FEM, and FHA are calorimeters, CMU, CMP, CMX, and FMU 
are muon chambers, BBC is a scintillation counter, and CPR and CES are multiwire proportional 
chambers. From Ref. g. 

penetrate the calorimeter. There are tradeoffs between these various systems: of the two de- 
tectors, CDF puts relatively more emphasis on tracking, while D0 emphasizes calorimetric 
measurements. 

We use a coordinate system centered on the detector, with the z-axis along the beam 
direction, and the x- and y-axes defining the transverse plane. We also use the polar angle 
9, the azimuthal angle 0, and the pseudorapidity r], defined as r\ = — ln(tan#/2). It is also 
common to calculate the angle between directions of two objects in terms of the distance 
between them in the (77, 0) plane, as R = J (A77) 2 + (A0) 2 . 



3.2.1. CDF 

At the core of the CDF detector is a cylindrical tracking volume contained within a 
large (4.8 m long by 1.5 m radius) superconducting solenoid that generates a field of 1.4 T. 
Immediately surrounding the beam pipe is the four-layer silicon vertex detector (SVX). This 
51 cm long device provides precise track reconstruction in the plane normal to the beam, 
measuring track impact parameters with a resolution of a = (13 + 40 (GeV/c)/pr) /mi. 
This is sufficiently precise to identify displaced vertices from b- and c-quark decays. The 
luminous region at CDF has a la width of about 30 cm, giving the SVX a geometrical 
acceptance of about 67%. The original SVX detector suffered significant radiation damage 
during Run la, but was replaced before the start of Run lb with a similar detector equipped 
with radiations electronics^ Outs.de the SVX is a vertex drift chamber (VTX), which 
is used to measure the z position of the interaction vertex to a precision of 1 mm. The VTX 
in turn is mounted inside of the central tracking chamber (CTC), a large cylindrical drift 
chamber which measures the curvature of tracks passing through the magnetic field. 

Surrounding the tracking volume are the calorimeters. There are three distinct calorime- 
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Fig. 10. The D0 detector, in both isometric and cross sectional views. CF, EF, and SAMUS are 
the three regions of the muon system (the small angle SAMUS system is not used in top analyses). 
From Ref. |l|. 

ter regions: central (\rj\ < 1.1), end-plug, and forward. In each region, there is an electro- 
magnetic calorimeter backed by a hadronic calorimeter. The electromagnetic calorimeters 
use lead as an absorber, while the hadronic calorimeters use iron. The active media are 
scintillator tiles in the central region and gas proportional chambers in the end-plug and 
forward regions. The tower geometry is projective — that is, the towers point to the nominal 
interaction point — with a cell size in the central region of Arj x A<p « 0.1 x 0.25. A layer 
of proportional wire chambers is located at shower maximum in the central electromagnetic 
calorimeter, and additional proportional chambers located between the solenoid and the 
CEM sample the early development of electromagnetic showers in the solenoid. The energy 
resolution, a(E)/E, of the calorimeters in the central region is about 13.7%/ WE © 2% for 
electromagnetic showers and about 50%/a/E © 3% for single hadrons. The thickness of the 
central calorimeters is about 18 radiation lengths for the electromagnetic section and 4.5 
interaction lengths for the hadronic section. 

Muons are identified using drift chambers surrounding the calorimeter, the central muon 
system and central muon upgrade. There is about 0.6 m of steel between these two sets of 
chambers. Together, they provide coverage out to about \r]\ < 0.6. Coverage is extended 
out to about \rj\ < 1.1 by the chambers of the central muon extension. 

3.2.2. D0 

D0 is the newer of the two detectors, commissioned during the summer of 1992. Its 
tracking volume is relatively compact (1.35 m long by 0.78 m radius) and nonmagnetic. 
Nested around the beam pipe are the vertex drift chamber (VTX), a transition radiation 
detector (TRD), and the central drift chamber (CDC). The tracking volume is capped on the 
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ends by two forward drift chambers (FDCs). The trajectories of tracks of charged particles 
can be measured with a resolution of 2.5 mrad in and 28 mrad in 9, and the ^-coordinate 
of the interaction vertex can be measured with a resolution of about 8 mm. The central 
tracking system also measures the ionization of tracks in order to distinguish between single 
charged particles and e + e~ pairs from photon conversions. 

The calorimeter is divided into three parts: the central calorimeter (CC) and the two 
end calorimeters (ECs). They each consist of an inner electromagnetic (EM) section, a fine 
hadronic (FH) section, and a coarse hadronic (CH) section. The absorber in the EM and 
FH sections is depleted uranium; in the CH section, it is a mixture of stainless steel and 
copper. The active medium in all cases is liquid argon. 

The EM sections of the calorimeters are about 21 radiation lengths deep, and are read 
out in four longitudinal segments (layers). The hadronic sections are 7-9 interaction lengths 
deep, with either four (CC) or five (EC) layers. The transverse segmentation is pseudo- 
projective (that is, although each cell is nonprojective, they form towers which are), with 
a cell size of Arj x A<p = 0.1 x 0.1. In the third layer of the EM calorimeter, near the 
shower maximum, the segmentation is twice as fine in each direction, with a cell size of 
At/ x A(f) = 0.05 x 0.05. The energy resolution is about 15%/a/E ©0.4% for electromagnetic 
showers and 50%/\/E for single hadrons. The resolution is substantially worse, however, in 
the transition regions between the CC and the ECs, due to the presence of a large amount 
of uninstrumented material. Some of the energy that would otherwise be lost is collected 
in extra argon gaps mounted on the ends of the calorimeter modules ( "massless gaps" ) and 
in scintillator tiles mounted between the CC and EC cryostats (intercryostat detectors, or 
ICDs). 

The D0 muon system consists of a "wide angle" system, covering |r/| < 1.7, and a 
"small angle" (forward) system, extending the coverage out to \r)\ < 3.3. For studying top 
quark decays, D0 uses only the wide angle system. This system consists of four planes of 
proportional drift tubes in front of magnetized iron toroids, with a magnetic field of 1.9 T, 
and two groups of three planes of proportional drift tubes behind the toroids. The magnetic 
field lines and the wires in the drift tubes are oriented transversely to the beam direction. 
The muon momentum is measured from the muon's deflection angle in the magnetic field of 
the toroid. The total amount of material in the calorimeter and iron toroids varies between 13 
and 19 interaction lengths, making the background from hadronic punchthrough negligible. 
In addition, the compact central tracking volume reduces backgrounds to prompt muons 
from in-flight decays of n and K mesons. During Run lb, the forward muon chambers 
suffered radiation damage that reduced their efficiency Midway through the run, however, 
the damage was repaired. As a result, the D0 top quark analyses do not use forward muons 
for the first half of Run lb. 



3.3. Particle Identification 



This section summarizes the algorithms used by the two experiments to identify the 
various final-state objects in candidate tt events. For more details, see Refs. p||52|, and ^3|. 
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3.3.1. Quarks and Gluons 

As a quark or gluon leaves the site of a hard scattering it cannot remain free, but instead 
hadronizes (or fragments) into a collection or jet of (colorless) hadronic particles. This 
collection tends to lie in a cone around the direction of motion of the original parton, and will 
show up in a calorimeter as an extended cluster of energy. In order to compare measurements 
with theoretical predictions it is necessary to have a precise definition of a jet: that is, one 
must specify how calorimetric energy depositions (cells) are to be clustered into jets. This 
algorithm is, in principle, arbitrary. However, at hadron colliders, it is conventional to define 
jets by taking all calorimeter cells which lie within a cone of fixed radius R in the (rj, 0) plane. 
This choice is convenient because jets are approximately circular in these variables; further, 
the i?-width of jets of a given is independent of the jet rapidity. More importantly, this 
definition can be readily implemented in phenomenological calculations, thereby facilitating 
the comparison of theory with experimental data.EI 

In principle, not only is the jet algorithm arbitrary but also the cone radius R. In prac- 
tice, the choice of cone radius involves several competing considerations. Jets are extended 
objects, composed of a collection of particles from hadronization of the progenitor parton. 
The jet will be further broadened as the particles undergo showering in a calorimeter. Con- 
sequently, if R is too small, a substantial portion of the energy from the progenitor parton 
will lie outside of the jet cone. This effect can be corrected for on average. However, the 
smaller the cone radius, the larger the energy correction that must be applied and, therefore, 
the worse the energy resolution of the corrected jet. On the other hand, if R is made too 
large, one cannot resolve the energy depositions arising from closely spaced partons; instead, 
the depositions get merged together into a single jet. This is of particular concern for it 
events, which tend to have many jets in the final state. The optimum choice for R for it 
physics depends somewhat on the structure of the calorimeter, but appears to be around 
0.4-0.5: CDF chooses R = 0.4, and D0 uses R = 0.5. 

Although both experiments have a resolution for single hadrons that scales as a(E)/E ~ 
50%/a/E, the resolution achievable for jets is typically a(E)/E ~ 100%/\AB. Most of the 
particles comprising a jet are of relatively low energy, in which region nonlinear effects in 
calorimeter response become important. Jet resolutions are also degraded by effects such as 
gluon radiation, differences in calorimeter response to hadrons and electrons, energy falling 
out of the jet cone, and contamination from hadrons from the underlying event. 

The measurement of jet energies is subject to numerous systematic effects, for which one 
must correct. These include: 

• The intrinsic response of the calorimeter to jets. 

• Calorimeter nonuniformities, and regions with uninstrumented material (such as cracks 
between modules). 

• Energy from the underlying event and, at D0, noise from the radioactive decay of the 
uranium absorber. 

• QCD radiation of gluons outside of the jet cone. 

• The spreading of particle showers outside of the jet cone in the calorimeter. 
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The procedures involved in performing these corrections are quite complicated;BE3 we 
shall therefore only summarize the strategies used. 

• Dijet events, in which the transverse energies of the two jets should balance, can be 
used to calibrate one region of the detector relative to another, better characterized, 
region. 

• Events with an electromagnetic cluster recoiling against a jet can be used to calibrate 
hadronic calorimeters relative to electromagnetic calorimeters. 

• The absolute scale of electromagnetic calorimeters can be determined by comparing 
electron energies to their momenta measured in the tracking system (at CDF) or by 
using the known masses of resonances such as the Z, the 7r°, and the J/ip (at D0).E3 

• Contributions from the underlying event and noise can be studied using Monte Carlo 
simulations, and by comparing data taken under differing trigger conditions and lumi- 
nosities. 

• The broadening of showers in calorimeters can be studied using test beams.il Monte 
Carlo simulations are used to model the distribution of particles produced during 
hadronization of partons. 

CDF and D0 apply jet corrections at different points in their analyses. This should 
be kept in mind when comparing selections involving jet energies. CDF does not use the 
jet corrections for measuring cross sections (except in the all-jets channel), but does apply 
them for the mass measurement, after the event sample has been selected. D0, on the 
other hand, applies most corrections before making any analysis selections. The corrections 
include effects of jets spreading in the calorimeter, but not of particles originating from 
gluons radiated outside of the jet cone. D0 applies an additional correction in its mass 
analysis to include this effect. 

It is important to realize that there is not necessarily any one-to-one correspondence 
between quarks and gluons in the final state of the hard scattering and the detected jets. 
A jet may have insufficient energy to be selected as a jet (a typical requirement is that the 
jet energy be at least 15 GeV), or two partons may be sufficiently close together that their 
energies are merged together during jet reconstruction. Conversely, if a parton radiates a 
gluon with a large relative transverse momentum, then that gluon may be reconstructed 
as a separate jet. Moreover, nonclassical effects, such as partonic interference, are always 
present and place a fundamental limit on the validity of identifying a given jet with a specific 
progenitor parton. 

3.3.2. Electrons 

Electron identification is based on finding isolated clusters of energy in the electromag- 
netic sections of the calorimeter, along with a matching track in the central detector from 
a charged particle. Additional requirements are then made to further suppress background 
from QCD jets. The exact requirements vary between experiments and among different 
analyses. However, typical requirements are: 
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• The fraction fg of the cluster energy in the electromagnetic sections of the calorimeter 
should be > 90%. 

• The shape of the cluster should be consistent with expectations from test beams or 
Monte Carlo simulations. 

• The track momentum p should be consistent with the cluster energy E, E/p ~ 1 (CDF 
only). 

• The distance between the cluster and the extrapolated track should be small. 

• A track ionization consistent with a singly charged particle (D0 only). Photon conver- 
sions into e + e~ pairs typically deposit twice the charge expected from one minimum 
ionizing particle. 

• Transition radiation information consistent with an electron (used in D0 dilepton 
analyses) . 

• Isolation, based either on calorimetry or (at CDF only) on tracking. One typically 
requires that the energy or momentum in the region around the electron candidate be 
small. 

CDF accepts electrons in the central calorimeter with \r]\ < 1.0. For the dilepton channels, 
one electron may also be in the plug calorimeter, thereby providing extra coverage in the 
range 1.20 < \r)\ < 1.35. D0 accepts electrons out to \r]\ < 2.5 for the dilepton channels and 
out to \r)\ < 2.0 for the lepton+jets channels (although the efficiency is poor in the transition 
region between the central and end cryostats, 1.1 < \rj\ < 1.5). 

3.3.3. Muons 

At CDF, muons are identified by requiring a match between a CTC track and a track 
segment in the muon chambers. This provides coverage out to \q\ < 1.0. For dilepton 
channels, CDF also identifies muons in sections of the detector where there is no coverage 
from muon chambers but good central tracking. This extends coverage out to \q\ < 1.2, 
and fills in coverage of azimuthal holes for < 1.0. Additional selections are made on the 
following variables: 

• The energy deposited in the calorimeters along the muon track. 

• The distance of closest approach of the muon track to the beam line. 

• The distance in z between the interaction vertex and the muon track. 

• The distance between the extrapolated CTC track and the track segment in the muon 
chambers. 



• Isolation, defined in a similar manner as for electrons except that calorimeter-based 
isolation is used when there is no matching track in the muon chambers. 
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D0 identifies muons by tracks in the muon chambers. A matching central detector track 
is not required, but if one is present, it is taken into account in the momentum measurement. 
The following additional criteria are used for D0's muon selection: 

• The line integral over the magnetic field > 2.0 T • m. This rejects muons which pass 
through the thin portion of the toroid around \rj\ ~ 0.9. Such muons bend less so their 
momenta are poorly measured. They also have a larger contamination from hadronic 
punchthrough. 

• The energy deposited along a muon track in the calorimeter must be at least that 
expected for a muon (~ 2-3 GeV). 

• The impact parameter between the muon track and the interaction vertex must be 
small. 

• Muons are defined as either isolated or nonisolated, depending on whether the distance 
in the rj — <ft plane between the muon and the closest jet is greater than or less than 
0.5, respectively. 

For the first portion of Run lb, when the forward muon chambers were affected by radiation 
damage, the muon selection was restricted to the central region, with \r]\ < 1.0. For other 
run periods, muon acceptance extended out to \r]\ < 1.7. 

3.3.4. t Leptons 

Tau leptons are difficult to identify. A r will decay into an electron or a muon about 
36% of the time, and the only observable difference between the decays W — ► Iv and 
W — ► tv — > ivvv is that the lepton spectrum is somewhat softer in the latter case. If the r 
decays hadronically, it can be detected as a narrow jet with either one or three associated 
tracks. The branching ratios for these "one-prong" and "three-prong" hadronic decays are 
about 50% and 14%, respectively! The challenge is to reject the large background from quark 
and gluon jets. CDF manages this with two complementary techniques, one "track-based" 
and the other "calorimeter-based. "@ 

The track-based technique is sensitive only to one-prong r decays. It starts by finding 
an isolated, high-py (> 15 GeV/c), central (\t]\ < 1.0) track. The isolation requirement is 
based on the sum of the px values of all tracks in a AR cone around the high-pr track. This 
discriminates between r leptons and jets. The energy in the calorimeter around the track 
is required to be consistent with the track momentum. In addition, candidates consistent 
with being electrons (large fractional energy in the electromagnetic calorimeter) or muons 
(energy deposition consistent with a minimum ionizing particle) are rejected. 

The calorimeter-based technique starts from a calorimeter cluster with \rj\ < 1.2 with 
either one or three isolated charged tracks pointing at it. (The tracks must have px > 
1 GeV/c and lie within a 10° cone around the centroid of the cluster.) About 73% of one- 
prong and 41% of three-prong decays are expected to contain 7r°s, which can be identified 
from their tt° — > 77 decays. The candidate is rejected if there are more than two 7r°s. The 
r lepton px is defined as the sum of the transverse momenta of all candidate tracks plus the 



27 



transverse momenta of all 7r°s; it must satisfy pt > 15 GeV/c. Also, the total invariant mass 
constructed from the tracks and the ir°s must be < 1.8 GeV/c 2 . The calorimeter cluster 
must be narrow, and its energy must be consistent with the total r momentum. Finally, 
clusters consistent with being electrons or muons are removed. 

The efficacy of the calorimeter-based selection is demonstrated in Fig. [□], which shows 
the track multiplicity for a monojet sample which required one jet with 15 < Ej- < 40 GeV, 
20 <]£ T < 40 GeV, and no other jet with Ex > 7 GeV. The excess in the one and three 
track bins from W — > rv is apparent; it is also seen that the calorimeter-based selection 
drastically reduces the background. About 45% of hadronic r decays satisfy the kinematic 
requirements; of these, about 55% are identified by this algorithm. 



3.3.5. Neutrinos 

Neutrinos do not interact in the detector with any significant probability, and therefore 
cannot be observed directly. Instead, their presence is inferred from an imbalance in the 
total transverse momentum in the event. Because the remnants of the spectator partons 
escape down the beam pipe, the longitudinal component of the neutrino momentum cannot 
be measured. 

The missing transverse energy $ T is defined by 

$t x = - Ei E i sin Qi cos 4>h (15) 

$T y = ~J2iEi sin 9i sin <f> it 

where the sums are over all calorimeter cells. If there are muons in the final state, then their 
transverse momenta should also be subtracted from $ T . (Some analyses make use of the $ T 
measured by the calorimeter without correcting for muons. This quantity will be denoted by 
$7^.) The resolution of $ T is usually parameterized in terms of the total transverse energy 
in the event (a scalar sum). CDF quotesi a resolution of o~{$ T ) = 0.7a/E Et (units in GeV), 
while D0 quotes!! a{$ T ) = 1.08 GeV + 0.019 ££ T . 



3.3.6. Tagging b-Jets 

A prominent feature of ti decays is that each event contains two 6-quarks. This is in 
contrast to the principal backgrounds, in which heavy flavors are expected to be relatively 
rare. Clearly, a method for identifying, or tagging, 6-jets would be valuable in separating 
tt events from background. The two strategies developed for doing this, soft-lepton tagging 
(SLT) and displaced- vertex tagging (SVX, using the silicon vertex detector), are discussed 
below. 

3.3.6.1. b- Tagging with Soft Leptons Approximately 22% of the time, a decaying b- 
quark will yield a muon, either directly or through a sequential decay via a c-quark. The 
branching ratio to electrons is also 22%. One method of 6-tagging is based on observing such 
leptons close to a jet. These leptons are much softer than the leptons from W boson decay, 
typically below about 20 GeV/c. (See Fig. |l2| .) Because they are soft and not isolated, 
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Fig. 11. Track multiplicity in the monojet data sample, a) No r identification requirements, 
b) All calorimeter-based r identification requirements, except for track multiplicity. From Ref. 58. 
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Fig. 12. spectra of leptons from decays of b- and c-quarks in it Monte Carlo events 
{m t = 160 GeV/c 2 ). From Ref. |. 
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their detection efficiencies are significantly lower than for the high-p^ isolated leptons from 
W boson decay. 

Both experiments tag 6-jets with this method. CDF uses both electron and muon tags, 
while, in the analyses reviewed here, D0 uses only muon tags. Without a central magnetic 
field, electron tagging is less effective at D0. However, even at CDF, the electron tags play 
a relatively minor role: the efficiency for electron tagging is about 1/3 of that for muon 
tagging. Another difference is that CDF selects soft leptons with p? > 2 GeV/c, while 
D0 requires px > 4 GeV/c because muons with lower momenta have insufficient energy 
to traverse the calorimeter and toroid. However, D0 has a larger muon acceptance than 
CDF. In addition, the fake rate for muon tags is significantly smaller at D0, due both to 
the large amount of material in the calorimeter and muon toroid (which reduces hadronic 
punchthrough) and the small size of the tracking volume (which reduces background due 
to in-flight tt/K decays). When all factors are accounted for, the probability of finding a 
lepton tag in a ti decay is nearly the same for both experiments, about 20%. 

3.3.6.2. b-Tagging with Displaced Vertices Another method of tagging 6-jets profits 
from the relatively long lifetime of b and c hadrons (about 1 ps). Given the typical boost of 
6-quarks in ti events, this allows the b hadrons to travel up to several mm before decaying. 
Detecting a vertex displaced by this distance is well within the capabilities of modern silicon 
microstrip detectors. 

At present, CDF is the only experiment to have operated a silicon vertex detector (the 
SVX) at a hadron collider.0 The tagging algorithm works by finding combinations of at 
least two tracks consistent with originating from a vertex displaced from the primary vertex 
of the hard interaction. Such displaced vertices are sometimes called secondary vertices. 
For each possible secondary vertex, one estimates the distance {L xy ) in the transverse plane 
between that vertex and the primary one, along with its associated uncertainty [phxy)- In 
order to be accepted as a &-tag, a vertex must satisfy the condition L xy /aL xy > 3. The sign 
of L xy is given by the sign of the dot product between the direction of L xy and the direction 
of the vector sum of the momenta of the tracks used. It is predominantly positive for real 
6-decays; displaced vertices with negative L xy are due primarily to track mismeasurements. 
Jets which contain many mismeasured tracks but no real secondary vertices are equally likely 
to have L xy positive or negative, a fact which is used to measure the background from this 
source. Only tags with positive L xy are accepted as 6-tags. The probability of finding at 
least one displaced vertex (or SVX) tag in a it event is (39±3)%H This includes a geometric 
efficiency of 67% caused by the length of the luminous region relative to the length of the 
SVX. See Fig. O for a sample result from this technique. 



3.4. Characteristics of Signal and Background Events 

We have discussed the decay modes of the top quark and the experimental signatures 
for it production in Sec. |2|. We have also outlined how the objects in the final state of a ti 
decay are identified and measured in the two detectors. Here, we shall discuss briefly the 
kinematic properties of ti events in the various decay modes, and how these properties differ 
for background processes. 

In the dilepton channels, the two leptons arise from W boson decays, and therefore 
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Fig. 13. The distribution of "pseudo decay length" A ~ ML xy /pT from CDF B — > J/ip — > /U/i 
data. The points are the data, and the curves are the signal and background models. The heavily 
shaded area is the background; it is asymmetric because of sequential semileptonic B decays in the 
dimuon sample. The lightly shaded area is the signal contribution, and the unshaded area is from 



prompt J/tfj mesons. From Ref. 6C. 
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Fig. 14. Expected distributions for e\i dilepton events of (a) lepton pr (two entries per event), 
(b) $ T , and (c) lepton rj (two entries per event). The solid histograms are it —* ep + X signal 
events (generated with herwig with m t = 175 GeV/c 2 for pp collisions at y/s = 1.8 TeV). The 
dashed histograms are Z + jets — ► tt + jets — > ep + jets events (also generated with herwig). All 
histograms are normalized to unity. We require that events have pj* > 10 GeV/c, $ T > 10 GeV, 
and at least two jets with Er > 15 GeV and \rj\ < 2.0. 



tend to be central in pseudorapidity, isolated, and of large transverse momenta. The two 
6-jets also have high transverse energies. There are two high-p^ neutrinos, so the missing 
transverse energy tends to be large as well. The major background for these channels is 
from the Drell-Yan process, which produces isolated lepton pairs in profusion. Additional 
jets can arise from initial or final state radiation. The Drell-Yan process yields ee and \ip, 
events directly, while ep events can be produced via rr production and subsequent decay. 
The ee and pp events at the Z resonance can be easily eliminated by rejecting events in 
which the invariant mass of the two leptons is consistent with that of the Z boson, Mz- 
Additional rejection of background is obtained by requiring at least two jets in the final 
state. In Drell-Yan events, the additional jets are due to gluon radiation; consequently, 
every additional jet reduces the cross section by a factor of 0(1/ a s ) ~ 6.7. Requiring a 
large ~$ T further reduces the background. In case of Z — > rr — > (lvv)(lvv), the Z boson 
cannot be reconstructed because of the presence of four unobserved neutrinos in the final 
state. However, the leptons in these events have much smaller than those in it events. 
So, requiring two high-p^ leptons, two or more jets, and large $ T will greatly suppress this 
background. Other backgrounds which must be considered are diboson production (W^W^, 
WZ, ZZ), QCD production of bb (with semileptonic decays of the 6-quarks), and QCD 
events with jets misidentified as leptons. The distributions of several kinematic quantities 
for dilepton events are shown in Figs. [14], [15], and |i~6| . 

The dominant background to the lepton+jets channel comes from the production of 
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Fig. 15. Expected distributions for e/i events of (a and b) the transverse energies of the two 
leading jets and (c) the jet r\ (two entries per event). See Fig. [l4| for further details. 
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Fig. 16. Expected distributions for e/i events of ffj- = X^jets + ^t- See Fig. 14 for further 
details. 
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Fig. 17. Expected distributions for lepton+jets events of (a) lepton px, (b) $ T , and (c) lepton r\. 
The solid histograms are it signal events (generated with herwig with m t = 175 GeV/c 2 for pp 
collisions at ^/i = 1.8 TeV). The dashed histograms are W+ > 4 jets events (generated with 
VECBOS). All histograms are normalized to unity. We require that events have p T > 15 GeV/c, 
$ T > 15 GeV, and at least four jets with E T > 15 GeV and \rj\ < 2.0. 



single W bosons in association with jets. Other backgrounds include false lepton events 
with mismeasured $ T and/or misidentified 6-tags, Wbb and Wcc production, and QCD 
heavy flavor processes. The distributions of E T , i]£, and $ T for the signal and for the 
VT+jets background are shown in Fig. [17|. The major differences between the signal and 
background processes stem from the number of jets in the event and the event kinematics. 
Ideally, a lepton+jets it event has four jets, two of which come from 6-quarks. However, 
as discussed in Sec. |3.3.1| , the number of observed jets can be greater or lesser than four. 
Therefore, the analyses usually require at least three jets if a 6-tagged jet is required in the 
event, and at least four jets otherwise. The jets in background events have lower transverse 
energies than those in the signal, and are produced over a wider range of rj, as is evident 
from Fig. [H| 

Certain variables describing the overall shape of the event provide powerful means of 
discriminating lepton+jets signal from background.i!0 One such variable is the aplanarityi! 
A, defined as 2/3 times the smallest eigenvalue of the total normalized momentum tensor 
in the event. This is 0.5 for spherical events and zero for planar or linear events. Top quark 
events are expected to be more spherical than radiative QCD processes, and hence to have 
larger aplanarities. Another useful variable is the sum of the transverse energies of all jets, 
Ht- This variable reflects the "temperature" of the interaction; a large Hx is a signature of 
the decay of massive objects.il Distributions of Ht and A are shown in Fig. [L9. 
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Fig. 18. Expected distributions for lepton+jets events of (a-d) the transverse energies of the 
four leading jets and (e) the jet r) (four entries per event). See Fig. [l?] for further details. 
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Fig. 19. Expected distributions for lepton+jets events of (a) Ht and (b) aplanarity. See Fig. 17 
for further details. 



4. DISCOVERY OF THE TOP QUARK 



The pursuit of the top quark began in earnest in the late 1970s, shortlvafter the discovery 
of the 6-quark at Fermilab.il Searches at the e + e~ colliders PETRaH (1979-84, y/s = 
12-46.8 GeV), TRISTAN§ (1986-90, yfe = 61.4 GeV), and SLCi and LEP§ (1989-90, 
y/s « Mz) eventually raised the lower limit on the top quark mass to 45.8 GeV/c 2 . (This 
limit and all others mentioned in this section are 95% confidence level results.) Owing to 
the pioneering work of S. van der Meer, C. Rubbia, and others, at CERN and elsewhere, 
high-energy pp colliders were developed in the 1980s. The first was the ISR (intersecting 
storage rings) at CERN§1. Next came the SppS, also at CERN. With yfs up to 630 GeV, 
this machine had a beam energy an order of magnitude higher than the ISR. This was 
followed by the Fermilab Tevatron, with y/s = 1.8 TeV. These machines provided much 
higher center-of-mass energies, enabling searches for particles with higher masses. 

Searches for the top quark at pp colliders do not provide direct limits on the mass of the 
top quark, but rather upper limits on its production cross section. These upper limits can 
be turned into lower limits on the mass usingcalculations of the production cross section. 

In 1984, the UAl collaboration reportedE2l evidence for the production of top quarks 
with m t = 40 ± 10 GeV/c 2 . In a subsequent analysis, however, with a larger data sample 
and a more thorough evaluation of backgrounds, the putative signal vanished!0 A limit of 
m t > 45 GeV/c 2 was inferred from this latter analysis. The UAl and UA2 experiments 
continued running through 1989, eventually setting limits of m t > 60 GeV/c 2 and m t > 
69 GeV/c 2 , respectively.^ Yet, even as they were in their last stretch, the CDF detector 
came online at the Fermilab Tevatron in 1988 and started recording pp collisions at the 
unprecedented yfs = 1.8 TeV, racing the CERN collaborations for evidence of top quark 
production with m t < Myy. The CDF collaboration soon set limits of m t > 77 GeV/c 2 
from the e + jets channel and m t > 72 GeV/c 2 from the e/i channel.^ Even with a smaller 
integrated luminosity (JCdt = 4.4 pb _1 for CDF vs. JCdt = 7.5 pb _1 for UA2), this limit 
was already better than could be achieved at the SppS experiments (because of the higher 
beam energy of the Tevatron). CDF later extended this analysis, adding the ee, /x/x, and 
/i + jets channels and using soft-lepton 6-tagging in the lepton+jets channels, arriving at 
a final limit from the 1988-89 runO of m t > 91 GeV/c 2 . Given these limits, the CERN 
experiments were out of the running; the search would continue only at the world's highest 
energy collider, the Fermilab Tevatron. 
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Collider operations at the Tevatron resumed in July, 1992, at which time the CDF 
detector was joined by the newly-commissioned D0 detector. Collider running continued 
through 1996, with the Tevatron reaching peak luminosities of over 2 x 10 31 cm~ 2 s _1 , a 
factor of 10 higher than the previous run, and twenty times the design luminosity The 
average luminosity for this period was 1.4 x 10 31 cm" 2 s _1 . 

Using the data from the first period of running (Run la), with JCdt = 13.5 pb -1 , 
the D0 experiment soon set a limit of m t > 131 GeV/c 2 using the ee, e/i, e + jets, and 
H + jets channels.il'0 (This limit was later revised downwards to 128 GeV/c 2 because of 
a recalibration of the luminosity at D0.) In April 1994, however, the CDF collaboration 
claimed!@ the first evidence for ti production. With an integrated luminosity from Run la 
of 19.3 pb -1 , CDF observed twelve candidate events in the dilepton and lepton+jets channels 
and estimated a 0.26% probability for the background to fluctuate to at least that many 
events. Assuming that the excess was due to ti production, the cross section was measured 
to be a t j = 13.9l 4 ; 8 pb. Under the same assumption, CDF also measured the top quark 
mass using the 6-tagged sample, obtaining a result of mt = 174± lOl} 3 , GeV/c 2 . Meanwhile, 
the D0 collaboration had reoptimized its analysis to search for high-mass top quarks, with 
m t ~ 180 GeV/c 2 . Nine candidate events were observed,il'@ compared to 3.8 ± 0.9 events 
expected from background. Taking m t = 180 GeV/c 2 yielded a it production cross section 
of 8.2 ± 5.1 pb. The chance of the observed signal being an upward fluctuation of the 
background was calculated to be 2.7%; therefore, D0 concluded that the excess was of 
insufficient statistical significance to demonstrate the existence of the top quark. 

Run la brought several spectacular it candidate events. In September, 1992, CDF 
recorded a beautiful it — > (W + b)(W~b) — > (e + ub)(qqb) candidate. A display of the SVX 
tracks in the event is shown in Fig. |20|. The event has an isolated electron, large $ T , and 
two jets with clearly identified displaced vertices (indicative of 6-quark decays). This event 
will surely find its way into the textbooks as an ideal top-antitop event! A kinematic fit of 
this single event to the ti decay hypothesis yields m t = 170 ± 10 GeV/c 2 . 

Another spectacular event is a dilepton (e/x) event recorded by D0 in January, 1993. 
An event display is shown in Fig. |2l[ This event has two high-p^ leptons (Ej. = 98.8 GeV, 
f T = 194.6 GeV/c), large $ T (100.7 GeV), and two jets (Eft = 26.1 GeV, Eft = 23.0 GeV). 
A multivariate Fisher discriminant analysis^! showed that this event is eighteen times more 
likely to be a it event than a [Z — > rr) + jets event, and ten times more likely to be ti than 
WW + jets. 

Run la continued through June, 1993. Over the summer, the accelerator was upgraded, 
with improvements to the linac and the installation of electrostatic separators. In Decem- 
ber, collider operations resumed with Run lb. By February, 1995, both experiments had 
quadrupled their data sets, and had observed large excesses of events over background that 
were fully consistent with the ti hypothesis. Finally, on March 2, 1995, the collaborations 
announced that the long search was over: the top quark had been found.SH 

CDF, in its 67 pb -1 data set, observed 37 lepton+jets events with at least one 6-tag. 
In this sample, there were 27 SVX 6-tags, compared to 6.7 ± 2.1 tags expected from back- 
ground, and 23 SLT 6-tags, compared to an expected background of 15.4 ± 2.0 tags. Six 
dilepton events were also seen, compared to 1.3 ±0.3 events expected from background. The 
probability for the estimated background to fluctuate to at least the observed number of 
events was calculated to be 1 x 10~ 6 , corresponding to a 4.8cx deviation for a Gaussian dis- 
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Fig. 20. A it —* e + jets candidate event from the CDF experiment. Jets 1 and 4 are identified 
by the SVX detector as 6-quark decays with vertices displaced from the primary event vertex. Jets 
2 and 3 could have come from a W boson, and the e + and v from the other W boson in the event. 
From Ref. fF?]. 
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Fig. 21. An e\i + jets candidate event from the D0 experiment. This is an end view of the 
tracking detectors and the central calorimeter. The muon is shown by the dotted line near the 
top, the electron is the large EM deposition near the bottom, and the two jets are the smaller 
depositions at about 7 and 8 o'clock. The $ T is shown by the narrow bar at about 5 o'clock. 
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tribution. CDF obtained a total tt cross section of 6.8l|; 4 pb, and a mass for the top quark 
of m t = 176 ± 8 (stat) ± 10 (syst) GeV/c 2 . 

Simultaneously, D0, using approximately 50 pb -1 of data, observed 17 events over an 
estimated background of 3.8 ±0.6 events. The probability for the background to fluctuate to 
at least the measured yield was 2 x 10~ 6 , equivalent to 4.6cx for a Gaussian distribution. D0 
measured a top quark mass of m t = 199^2? (stat) ± 22 (syst) GeV/c 2 and a it production 
cross section of 6.4 ±2.2 pb. 



5. MEASUREMENT OF THE tt PRODUCTION CROSS SECTION 

5.1. General Strategy 

Measuring the pp — > tt production cross section can be done in several decay modes. 
These are categorized as either dilepton, lepton±jets, or all-jets, as discussed in Sec. |2|. 
These can be further subdivided based on the lepton flavors in the final states, on whether 
or not 6-tags are present, and on the method used for 6-tagging. For each subchannel, the 
cross section is given by 

(N — B) 

°« = 7T-jd- (16 » 

where N is the number of observed events, B the estimated background count, e the total 
detection efficiency for it events in the subchannel, and JCdt the integrated luminosity of the 
data set. The total efficiency e includes the branching fraction, the geometrical acceptance 
of the detector, and the efficiencies for trigger selection, identification of leptons and jets, 
and kinematic selections. For any given set of kinematic selections, the signal efficiency 
depends on the mass of the top quark; therefore, the measured cross section depends on 
the assumed value of the mass. For some channels, the background prediction is based on 
data samples in which there is a small contribution from tt production. In such analyses, 
the measured tt cross section is used to correct the background prediction iteratively, until 
the measured cross section is stable. To obtain a final measurement of the cross section, all 
channels are combined, taking into account any correlated uncertainties, such as those on 
integrated luminosity and on the signal and background models. 

Both collaborations use the HERWIG0 Monte Carlo program as the primary model for tt 
events. In addition, CDF uses PYTHIA0 to assess systematic uncertainties due to modeling, 
while D0 uses isajetH The W / ±jets background, which is the principal background in the 
lepton±jets channels, is modeled using a combination of collider data and simulations based 
on the VECBOS programS The other important background in these channels arises from 
misidentification of one or more jets as leptons in multijet events. The misidentification rate 
(called the fake-lepton probability) is obtained from the data, as are the lepton identification 
efficiencies and the fe-tag probabilities of the background. Other (smaller) backgrounds 
are estimated using a combination of Monte Carlo simulations and object identification 
efficiencies measured from data. 
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5.2. CDF Analyses 

The most up-to-date analyses from CDF use data from Run la (1992-93) and Run lb 
( 1994-95). SliBH Since Run lc was short, CDF chose not to trigger on top-like events, but 
instead used the run to pursue other specialized studies. The total integrated luminosity 
used in the top analyses is 109 ± 7 pb -1 . If one assumes the predicted it cross section of 
5.5 pb (at m t = 175 GeV/c 2 ), about 600 it events should be present in the CDF data. A 
subset of these data (JCdt ~ 67 pb -1 ) supported CDF's observation of the top quark.El 

Using the full data set, CDF has updated analyses in the dilepton and lepton+jets/b-tag 
channels and has performed new analyses in the r dilepton (er and fir) and all-jets (all 
hadronic) channels. The r dilepton channels are of special interest because, if charged Higgs 
bosons H ± with m H < mt exist, they could produce an excess of events in these modes via 
the decay chain t — > H + b — > r + u T b. 

The all-jets final state accounts for 44% of all it events. It is therefore both important 
and interesting to test this key prediction through an independent observation of signal in 
this channel. 

CDF calculates o~ t j assuming m t = 175 GeV/c 2 , while D0 uses m t = 172 GeV/c 2 . 
This corresponds to a typical difference of ~ 0.3 pb in extracted cross sections, due to the 
dependence of the efficiency on mass. 



5.2.1. Dilepton Channels 

We discuss first the analyses in the "standard" dilepton channels, ee, e/i, and and 
leave dilepton channels with identified r leptons for later. 

The initial event selection in the standard dilepton analyses requires the presence of 
two oppositely charged high-p-r leptons (e or fi), two or more jets as expected from the 
6-quarks, and large $ T as the signature for the neutrinos. The kinematic selection criteria 
are shown in Table Since both leptons in dilepton it events come from W boson decays, 
at least one of the two leptons is required to be isolated. The criterion for lepton isolation 
is that the transverse energy in the calorimeter in a cone of R = 0.4 around the lepton 
be less than 10% of the lepton's E T (or p T ). The dominant background to the ee and \x\x 
channels comes from Z(— > it) + jets. This is largely eliminated by rejecting events with 
a dilepton invariant mass M# within a narrow window about the Z boson mass, that is, 
with \Mu — Mz\ < 15 GeV/c 2 . Events containing an isolated photon with Ex > 10 GeV 
consistent with a radiative decay of a Z boson are removed. The requirement on the $ T is 
tightened (to $ T > 50 GeV) when the Ift T vector is nearly collinear with either a lepton or 
a jet (A(j)(]£ T ,£ or jet) < 20°). This suppresses backgrounds from Z — > rr, where the two 
r leptons (and hence their decay products) are spatially close when the decaying Z boson 
has high momentum, and background from the Drell-Yan continuum where $ T arises from 
mismeasurements of jet or lepton energies. The latter process has a very large cross section, 
and therefore is an important source of background. The distribution of A0(^ T , i or jet) 



vs. $ T is shown in Fig. [22| for events that pass all but the final selection criterion. The 
expected distribution for it events, calculated with herwig assuming m t = 175 GeV/c 2 , is 
superimposed for comparison. Nine events — seven e/i, one ee, and one \x\x — survive all 
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Table 5. Basic kinematic requirements in various channels in the CDF cross section measure- 
ments. 





Standard Dilepton 


r-Dilepton 


I + jets/6-tag 


All jets 




(ee, eu, uu) 


(er, lit) 






Lepton px{GeV / c) 


> 20 


V T rV > 15 

p e J" > 20 


> 20 




Lepton \rj\ 


< 1.0 


\r] T \ < 1.2 
|r/ e ^| < 1.0 


< 1.0 




$ T (GeV) 


> 25 




> 20 




S $t (GeVV) 




> 3 






Jet E T (GeV) 


> 10 


> 10 


> 15 


> 15 


Jet \q\ 


< 2.0 


< 2.0 


< 2.0 


< 2.0 


Number of jets 


> 2 


> 2 


> 3 


> 5 


H T (GeV) 








> 300 


Hf{GeV) 




> 180 







the requirements. The distribution of H^ 1 (= Vt + Pt + $t + &t + ) f° r candidate 



events is compared to expectations for signal and background in Fig. |23 



The estimates of the backgrounds from various sources are listed in Table || After the 
initial event selection, Drell-Yan (Z/j* — > ee,///i) production continues to be the major 
background for the ee and channels. The e/x channel, the cleanest of the dilepton chan- 
nels, has background mainly from decays of W boson pairs (PFPF) and Z — > tt. In all cases, 
additional jets can arise from gluon radiation, and ]£ T from either neutrinos or mismeasure- 
ment of energies. CDF also estimates the background from false signatures due to particle 
misidentification, such as a jet or a track faking one of the leptons, or overestimated Ift T due 
to mismeasured muon momenta. The backgrounds from Drell-Yan production, fake leptons, 
and mismeasured tracks are estimated from the data. Other backgrounds are calculated 
using Monte Carlo simulations, which use lepton identification efficiencies estimated from 
data. 

For the acceptance of tt events, CDF takes the average of the results from herwig and 
pythia. The lepton identification efficiencies are measured to be 91% for muons and 83% 
for electrons, using Z — > £ + £~ data. (These efficiencies do not include the geometric accep- 
tance of the detector or the isolation cuts.) The overall efficiency for detection of tt events 
(including the branching ratio) is estimated to be e = (0.74 ± 0.08)% for m t = 175 GeV/c 2 . 
The uncertainty in the efficiency reflects uncertainties in event modeling (estimated from the 
differences between herwig and pythia) and in the simulation of the detector. Assuming 
a tt production cross section of 5.5 pb, CDF expects to see a total of 4.4 tt events in the 
standard dilepton channels, of which (58 ± 2)% are expected to be e/i, (27 ± 1)% /j/j, and 
(15 ± 1)% ee. Both charged leptons directly come from W boson decays in (86 ± 2)% of 
the observed dilepton it events. In the remaining events, one of the leptons comes from the 
decay chain W — > r — > ejii. The total tt acceptance increases by 35% as m t increases from 
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Fig. 22. The azimuthal angle A</> between the $ T vector and the nearest lepton or jet vs. $ T 
in CDF dilepton events. The small points show the distribution expected from it signal with 
rat = 175 GeV/c 2 . The larger symbols represent data. The dashed line shows the cuts on $ T . 
From Ref. 
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Fig. 23. Comparison of for the CDF dilepton candidate events (histogram) with the ex- 
pectation for tt production (m t = 175 GeV/c 2 ) plus background (lighter shaded area). The darker 
shaded area represents background alone. From Ref. pU 



Table 6. Expected backgrounds to the standard dilepton channels in the CDF data, corre- 
sponding to JCdt = 109 ± 7 pb- 1 . From Ref. ||. 



Background type 




Expected Number of events 
(All dilepton channels) 


Drell-Yan 




0.61 ±0.30 


Z — > TT 




0.59 ±0.14 


Fake leptons 




0.37 ±0.23 


WW 




0.36 ±0.11 


Mismeasured muon 


tracks 


0.30 ± 0.30 


QCD bb 




0.05 ±0.03 


Other (radiative Z, 


Wb, WZ, ZZ) 


0.10 ±0.10 


Total 




2.40 ±0.50 
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150 to 200 GeV/c 2 . 

With nine dilepton events in the data sample, an estimated background of 2.4 ± 0.5 
events, and an overall it efficiency of (0.74 ± 0.08)%, the cross section for it production for 
m t = 175 GeV/c 2 is found to be 8.2l 3 ; 4 pb. Of the nine candidate events, four have one 
jet tagged as a 6-jet by the SVX algorithm; of these four, two are also tagged by the SLT 
method. If all nine candidates are assumed to be from background processes, the number 
of expected 6-tagged jets is only 0.7 ± 0.2. 

CDF also searches for it events in the er and fir decay channels.il The total branching 
ratio for these channels is 4/81, the same as for all the standard dilepton modes (see Sec. f2.2| ). 
This could, in principle, double the efficiency for dilepton modes. That is not the case, 
however, since only hadronic decays of the r are used (branching fraction rs 64%: 50% one- 
prong and 14% three-prong decays), and the identification of r leptons is far less efficient 
than of electrons or muons. 

Recall from Sec. 3.3.4 that CDF uses two different r identification algorithms, one "track- 
based" and the other "calorimeter-based." Separate cross section analyses are performed for 
each of these selections, starting from samples containing events with r candidates along with 
oppositely charged leptons satisfying the kinematic cuts summarized in Table |^. Each event 
must also contain at least two jets (since we expect two 6-quarks) with Et > 10 GeV and 
\r)\ < 2.0, and have significant $ T due to the unobserved neutrinos. Instead of directly ap- 
plying a cut on $ Tl CDF applies a cut on the $ T significance, defined as = $ T / y/J2 E T 



for er events, and = $ T /\J J2 Et + Vt f° r ^ T events. Here J2 Et is the scalar sum 

of the transverse energies measured in the calorimeter and provides a measure of the Ifi T 
resolution. The requirement on the $ T significance is S-^ > 3 GeV 1 / 2 . The distributions of 

Sjp vs. ifjrp for the er and fir data samples as well as for it Monte Carlo events are shown 

in Fig. H A cut Hf > 180 GeV is also imposed (where Hf = p T + p T T + $ T + Ejets 
i = e or //). Using either r identification method, the same four events are left after all cuts, 
two er and two [it (shown as stars in Fig. ^4]). Three of these events are 6-tagged, and one 
has an SLT-SLT double-tag. 

The r identification and signal efficiencies are estimated using pythia it events decayed 
with the tauola package,ii which properly treats the r polarization. The r identification 
is cross checked using a data sample enriched in W — > tv t decays (see Fig. [11]). Of all the 
selected one-prong events, 19% and 38% are found only by the track-based and calorimeter- 
based techniques, respectively, and 43% by both. The overall signal efficiency is estimated 
to be e = (0.085 ± 0.010 (stat) ± 0.012 (syst))% for the track-based analysis, and e = 
(0.134 ± 0.013 (stat) ± 0.019 (syst))% for the calorimeter-based analysis. 

The systematic uncertainty is dominated by uncertainties in identification efficiencies for 
the r (6%) and e//i (7%) and in the hadronic energy scale of the calorimeter (5%). The 
uncertainty in the top quark mass contributes another 6%. If one uses the it production 
cross section of a t j = 7. pb, measured by CDF from other decay channels, one expects 
0.7 ± 0.2 (stat) ± 0.1 (syst) and 1.1 ± 0.3 (stat) ± 0.2 (syst) it events in the track-based 
and calorimeter-based analyses, respectively. So, if the top quark decays as predicted by the 
Standard Model, these channels are not expected to improve significantly the acceptance for 
signal. 
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Fig. 24. The distribution of S^. vs. $y for CDF events with a primary lepton and a r candi- 
date. Three of the four final candidate events (stars) have 6-tagged jets. From Ref. 
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Table 7. Numbers of events expected and observed for the CDF r-dilepton analyses. The 
one-jet columns are background dominated; the columns with two or more jets correspond to the 
signal region. From Ref. 



Selection 






Track- 


based 




Calorimeter-based 


ATjet > 10 GeV 




1 




> 2 




1 


> 2 


t fakes 


0.14 


± 


0.01 


0.25 ± 0.02 


0.47 


±0.03 


0.78 ±0.04 


Z/'y — ► t + t~ 


0.22 


± 


0.12 


0.89 ±0.28 


0.54 


±0.16 


1.48 ±0.38 


WW, wz 


0.14 


± 


0.06 


0.14 ±0.08 


0.20 


±0.09 


0.24 ±0.10 


Total Background 


0.50 


± 


0.14 


1.28 ±0.29 


1.21 


±0.28 


2.50 ±0.43 


Expected from it 


0.08 


± 


0.02 


0.70 ± 0.30 


0.13 


±0.03 


1.10 ±0.40 


Observed events 




1 




4 







4 



The primary background in the er and fir decay channels comes from Z/j* — > t + t~ ±jets 
events, where one r decays leptonically and the other hadronically. The other backgrounds 
include WW and WZ production and "fake r's" in W + jets events, where a jet is misiden- 
tified as a r. The fake r background is calculated by weighting the E? spectrum of all 
jets that could be misidentified as r's in the W+ > 3 jets sample by the fake rate, de- 
termined by applying the r selection criteria to a multijet event sample. The background 
estimates are shown in Table [7|. The total background for events with Nj et > 2 is esti- 
mated to be 1.28 ± 0.29 events for the track-based selection and 2.50 ± 0.43 events for the 
calorimeter-based selection. The measured it production cross sections, based on the four 
events observed, are a t j = 29. ll^ (stat) ±4.7 (syst) pb for the track-based selection and 
a t j = 10.2^\q^ (stat) ± 1.6 (syst) pb for the calorimeter-based selection. Unfortunately, 
because of their large uncertainties, these measurements do not improve the precision of 
the overall cross section measurement. But the analysis is important because non-standard 
decays of top, specifically t — > H + b with H + — > r + z/ r , could have produced an excess above 
the Standard Model prediction. No significant excess is evident. 



5.2.2. Lepton+Jets Channels 

To measure the it production cross section in the lepton±jets channels ,H CDF accumu- 
lated data using inclusive electron and muon triggers that required y\ > 18 GeV/c. A $ T 
trigger was used to compensate for small inefficiencies in the lepton triggers. 

The signal in this channel, it — > (W + b)(W~b) — > (£vb)(qq'b), is subject to a large 
background from VF±jets production. Since it is hard to identify W bosons that decay 
hadronically, these it events comprise a subset of the W(— > Iv) + jets sample. This sample 
is used for systematic studies of the background involving a W boson, as well as for extracting 
the signal. 

The initial data set consists of high-pr inclusive lepton events that contain a central 
< 1.0) isolated electron or muon with px > 20 GeV/c. The lepton isolation / is defined 
as the ratio of the transverse energy in a cone of R = 0.4 around the lepton (excluding the 
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lepton energy) to the pt of the lepton; the analysis requires I < 0.1. Events consistent with 
Z — > ££ or Z — > il^ are rejected, as in the dilepton analysis. From this inclusive lepton 
sample, events with $ T > 20 GeV are selected for further study. Events that pass the 
dilepton selection criteria are removed in order to keep the two samples disjoint. This gives 
an inclusive sample of W + jets events. 

Nominally, one expects to observe four jets per tt event in the lepton+jets channels. But, 
as discussed in Sec. |3.3.1| , the number of jets may be larger or smaller than this. Therefore, 
the analysis requires a minimum of only three jets. The kinematic selection criteria are 
summarized in Table [5]. 

CDF uses both the SVX and SLT algorithms for 6-tagging, with efficiencies of (39 ±3)% 
and (18 ± 2)%, respectively, for tagging at least one 6-quark in a tt event with at least three 
jets. These tagging efficiencies are estimated from a tt Monte Carlo simulation. After all 
cuts are applied, the W+ > 3 jets sample contains 34 SVX-tagged events containing a total 
of 42 SVX-tagged jets, and 40 SLT-tagged events with 44 SLT-tagged jets. Of these events, 
11 are 6-tagged by both algorithms. 

The primary background is from PF+heavy flavor production processes. For the SVX- 
tagged sample, CDF estimates this background by using vecbos and herwig to predict, 
as a function of jet multiplicity, the fraction of W + jets events that contain heavy flavor. 
This is then combined with the tag probability and applied to the observed W + jets events. 
Backgrounds from processes that do not produce a W boson are calculated from the data 
by measuring the tagging rate in a sample with low lfi T and without isolated leptons, where 
there are few W bosons. This is then used to predict the amount of contamination in the 
signal region. The background due to mistags can be found by assuming that the distribution 
of decay lengths from this source is symmetric around zero. CDF measures the negative half 
of this distribution from jet data and uses that to predict the number of events fluctuating 
above the cut. The background expected from single top production is also calculated, using 
the latest theoretical cross section® and efficiencies derived from Monte Carlo simulations: 
Additional small backgrounds from WW, WZ, and Z — > tt processes are calculated from 
Monte Carlo simulations. 

For the SLT-tagged sample, the background is mainly from W + jets events with fake 
leptons, K decays in flight, electrons from photon conversions, and events containing heavy 
flavor jets (Wbb, Wcc). These backgrounds are calculated from jet data by measuring the 
tag probability per track as a function of the track py. These tag probabilities are then 
applied to the tracks in the W + jets events to estimate the background from these sources. 

The signal efficiency, which includes the trigger, tagging and lepton identification ef- 
ficiencies, kinematic and geometric acceptances, and the branching fraction, is estimated 
using a combination of data and Monte Carlo events. For writ = 175 GeV/c 2 , the overall 
efficiency is (3.7 ± 0.3)% with SVX-tagging and (1.7 ± 0.3)% with SLT-tagging. 

The main systematic uncertainties arise from the jet energy scale (±5%), modeling of 
initial and final state gluon radiation (±2% and ±5%, respectively), Monte Carlo generators 
(±5%), detector resolution effects (±2%), and the instantaneous luminosity (±1%). 

The results of the analysis are summarized in Table ||; the SVX-tagged results are also 
displayed in Fig. |25|. Note that the SVX-based analysis achieves a signal to background ratio 
about five times better than the SLT-based analysis, demonstrating the clear advantage of 
6-tagging with a silicon vertex detector. Finally, for m t = 175 GeV/c 2 , CDF obtains a cross 
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Fig. 25. The CDF jet multiplicity distribution of SVX-tagged lepton+jets events. Note the 
clear excess above background for events with three or more jets. From Ref. 
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Table 8. Summary of CDF acceptance factors, event yields, backgrounds, and measured cross 
sections for each analysis channel, for nit = 175 GeV/c 2 . 





Dilepton 


Lepton+Jets 


All-Jets 






SVX 


SLT 


SVX 


2 SVX 


Ctotal 


0.0074 ± 0.0008 


0.037 ± 0.005 


0.017 ±0.0003 


0.044 ± 0.010 


0.030 ± 0.010 


Obs. Events 


9 


34 


40 


187 


157 


Background 


2.4 ±0.5 


9.2 ±1.5 


22.6 ±2.8 


142 ± 12 


120 ± 18 


Expected it 


4.4 


31.0 


14.3 






°it (P b ) 


o n+4.4 


6-2i?; 7 


9 2+ 4 - 3 


9 fi+ 4 - 4 


11 5+ 7 - 7 



section of 6.2±i; 7 pb using the SVX analysis, and 9.2±f;| pb using the SLT analysis. 



5.2.3. All-Jets Channel 

The all-jets channel has the advantage of having a large branching fraction (« 44%) 
and a fully reconstructible final state (since no high-p^ neutrinos are present). However, it 
is subject to an overwhelming background from QCD multijet production, with about one 
thousand times the rate of tt production. But, by requiring large transverse energy in the 
final state as well as the presence of 6-quarks, the multijet background can be reduced to 
manageable levels. 

CDF recently reportedil the first observation of tt decays in the all-jets channel (referred 
to as the all hadronic channel by CDF). The collaboration performs two analyses, which 
require either a single or a double 6-tag, and, in both cases, finds an excess relative to the 
estimated background. 

The data sample was collected using a special multijet online trigger, specifically devel- 
oped for tt decays into the all-jets final state. This trigger requires four or more clusters of 
contiguous towers in the calorimeter, each with > 15 GeV and a total transverse energy 
of £ > 125 GeV. Offline, jets are reconstructed using a cone algorithm with R = 0.4, 



and jet energy corrections are applied (see Sec. |3.3.1| ). Selecting events with five or more 
jets and requiring J2 E T > 300 GeV yields 21,840 events. Events with high-p T electrons and 
muons (which would appear in the dilepton and lepton±jets channels) are removed. After 
this selection, the QCD background still dominates tt production by a factor of about 110. 
But when one demands at least one SVX 6-tag, the signal to background ratio increases 
to S/B Pd 1/20, in a sample of 1596 events. Thereafter, two approaches (referred to as 
Techniques I and II) are used to further reduce background. 

In Technique I, two new variables are introduced that together prove decisive in achieving 
the final S/B ratio of « 1/4: 

• The centrality C = Ht/V§, where y/§ is the invariant mass of the multijet system. 

• The aplanarity A, computed from the jet momenta. 
The following cuts are then made on these variables: 
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• C > 0.75. 

• A > -0.0025# T3 + 0.54; (H T3 = H T - E^ 1 - E$ t2 ). 

The motivation for the first cut is the observation that the jets from it decay are more central 
than those from QCD production (see Sec. [J). The D0 collaboration uses a similar variable, 
except that the Ht is divided by the sum of the jet energies (see Sec. |5T3| ). The second 
cut provides a more optimal way to separate signal from background than that afforded by 
independent cuts on the aplanarity and Ht- (Cuts of this kind in the (.4, i?;r)-plane were 
first used by D0.@H) 

After all cuts, 187 events containing 222 6-tags remain. This is a statistically significant 
excess over the estimated background of 164.8 ± 10.8 6-tags. The background is largely from 
QCD production of heavy quark pairs and from fake fe-tags. The probability for finding 
a 6-tag in a multijet sample, using the SVX algorithm, is parameterized in terms of the 
jet Et, T) and track multiplicity, and the event aplanarity. To estimate the background, 
each multijet event is weighted by its tag probability. The sum of these weights, for the 
multijet events that pass all cuts, gives the background estimate. But, because the multijet 
sample contains it events, the background estimate must be corrected to account for the it 
contribution in the multijet sample. This procedure is done iteratively. The tag probabilities 
are cross checked by comparing predicted 6-tagged distributions with those observed, using 
an independent multijet sample. 

The uncertainty of the background estimate is dominated by systematics arising from 
the dependence of the tag probability on the run conditions (such as instantaneous lumi- 
nosity), the event shape, and kinematics. In Table |8] we summarize the event yields and 
the background estimates for N^ ets > 5. CDF estimates the probability of the background 
fluctuating to at least the number of observed 6-tags as V = 1.5 x 10~ 3 , corresponding to 
three standard deviations for a Gaussian distribution. 

From the number of tagged events and the background estimate, corrected for the it 
content, CDF extracts the number of it candidates to be 10.4 ± 6.0 for Run la data and 
34.7 ± 16.1 for Run lb data. The total signal efficiency is estimated to be (9.9 ± 1.6)% for 
m t — 175 GeV/c 2 . The uncertainty in the efficiency is largely due to systematic uncertainties 
in the jet energy scale (9%), in modeling of fragmentation (9%) and gluon radiation (11%). 
The 6-tagging efficiencies for the two data periods are (38 ±11)% and (46±5)%, respectively. 
The it cross section is measured to be a t j = 9.6 ± 2.9 (stat)!^ (syst) pb. 

In Technique II, two or more 6-tags are required in each event. Requiring the second 
6-tag significantly reduces the acceptance. Most of this is recovered, however, by removing 
the centrality and aplanarity requirements. CDF observes 157 such events with > 5 jets, 
while background from QCD heavy flavor production and fake double tags is predicted to 
be 122.7 ± 13.4 events. The number of signal events is extracted from a likelihood fit of a 
sum of contributions from fake double tags, QCD heavy flavor background, and it signal, to 
the observed number of events as a function of jet multiplicity. The number of it candidates 
is found to be 5.9 ± 3.9 for Run la and 31.6 ± 16.4 for Run lb. The corresponding numbers 
of background events are 21.1 ±4.5 and 98.4 ± 17.3. The detection efficiency before tagging 
for it events is (26.3 ± 4.5)% for m t = 175 GeV/c 2 . The estimated efficiencies for tagging 
> 2 heavy flavor jets are (7 ± 6)% and (12 ± 2)% for the two data-taking periods. The 
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cross section is a t j = 11.5 ± 5.0 (stat)!^ (syst) pb. The probability for the background to 
produce at least the observed excess is V = 2.5 x 10~ 2 , which corresponds to two standard 
deviations for a Gaussian distribution. 

CDF combines the results from the two techniques to calculate an overall cross section, 
taking into account the correlations in the efficiencies for the two techniques as well as 
the overlap between the two data samples (for which 34 events are in common). The 
correlation coefficient between the two techniques is estimated to be p = 0.34 ± 0.13. For 
m t = 175 GeV/c 2 , the combined measurement is a t j = 10.1 ± 1.9 (stat)!^ (syst) pb. 

5.3. D0 Analyses 

The D0 collaboration has used its full Run 1 data set of 125 ± 7 pb -1 to measure the 
tt production cross section.!!! Here, we review the D0 measurements for the various decay 
channels. These cross sections are calculated for m t = 172.1 GeV/c 2 .0 

5.3.1. Dilepton Channels 

D0 used the standard dilepton channels ee, e/x, and /i/i to support the discovery in 
1995. Since then, the collaboration has added the inclusive ev + jets channel. The dilepton 
channels are characterized by two high-p^ leptons, two high-.&r jets, and substantial missing 
transverse energy. The selection criteria are summarized in Table ^|. Note the use of the 
variable H T = Ht + E T , where E T is the transverse energy of the leading electron. This 
variable is found to provide better discrimination between signal and background than Ht 
alone. After additional cuts (discussed below), designed specifically to suppress the back- 
ground from Z + jets events, the sample contains three e/i events, one ee event, and one /z/i 
event. The distributions of the Ht variable for signal, background, and candidate events 
are shown in Fig. ^ (before the cut on Ht). 

D0 sought to increase the it acceptance by using the u eu" channel to recover some events 
that fail the standard kinematic selection. This channel requires only a single electron, but 
has very high lfi T and jet cuts. It is populated by events from the channels ee, e/i, and 
e + jets, and yields four events that pass the selection criteria listed in Table |9[ We note 
that the significance of these data is comparable to that of the dilepton data; therefore, the 
ev channel makes a useful contribution to the cross section measurement. 

The main backgrounds in the dilepton channels are from Z boson, Drell-Yan, and dibo- 
son processes. In the ee channel, the large background from Z — > ee is reduced by tightening 
the Iftrp cut from 25 GeV to 40 GeV if the dielectron mass lies within 12 GeV of the mass of 
the Z boson. For the /x/i channel, the Z + jets background is also dominant. Because the 
Hji and ee channels have identical kinematics one could, in principle, apply the same cuts 
in both cases. However, D0 does not measure the momenta of muons as well as those of 
electrons. Consequently, the missing transverse energy is more susceptible to mismeasure- 
ments of muon momenta, and is therefore less useful as a variable to discriminate signal 
from background. Instead, events are fit to the Z{— > nji) + jets hypothesis; events with 
acceptable fits (Prob(x 2 ) > 1%) are rejected. The e\i channel, by contrast, is particularly 
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Fig. 26. Ht distributions of dilepton events (for mj = 170 GeV/c 2 ) in the D0 analyses. The 
location of observed events is shown by dark histograms. 



Table 9. Kinematic event selection criteria for various tt decay channels applied in the cross 
section analysis by D0. All energies are in GeV, rj is the pseudorapidity, Ht = J2 -^t* w ^h. 
> 15 GeV, and H T = H T + E T . 





Dilepton 


I + jets 


I + jets/fi 


eu 


Lepton pt 


> 15 


> 20 


> 20 


> 20 




> 20 (ee) 








Electron \r]\ 


< 2.5 


< 2.0 


< 2.0 


< 1.1 


Muon \r]\ 


< 1.7 


< 1.7 


< 1.7 




$T 


> 20 (en) 


> 25 (e) 


> 20 


> 50 




> 25 (ee) 


> 20 (fi) 






Jet E T 


> 20 


> 15 


> 20 


> 30 


Jet \rj\ 


< 2.5 


< 2.0 


< 2.0 


< 2.0 


Number of jets 


> 2 


> 4 


> 3 


> 2 


H!t 


> 120 (ee,e/x) 










> 100 (nn) 


> 180 


> 110 




A 




> 0.065 


> 0.040 




H/rp 




> 60 






\v 




< 2.0 






Tag muon 




veto 


Pt > 4 










AK([i,}et) < 0.5 




M|f 








> 115 
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clean: with the selection criteria given in Table || D0 obtains a signal to background ratio 
of 10 : 1. Half the background comes from Z — > tt events. 

The W(— > ev) + jets background, the main background in the ev channel, is very effec- 
tively reduced by the application of a large transverse mass requirement: M T V > 115 GeV, 

where M T V = \jl\p T \ \t>t\ ~ ^Pt ' Pt- The size of the background is estimated from data: one 
scales down the number of W+ > 2 jets events by the efficiency of the transverse mass cut, 
the latter determined from a simulation of W+ > 2 jets events using the VECBOS program. 

The systematic uncertainties mainly arise from the uncertain knowledge of the jet energy 
scale and Monte Carlo modeling, and are estimated to be between 10 and 15%. The 
channel has an additional uncertainty of 10%, due to the Z — > fifi kinematic fit. From these 
channels, D0 measures o~ t j{m t = 172.1 GeV/c 2 ) = 6.4 ± 3.3 pb. 

5.3.2. Lepton+Jets Channels 

The lepton + jets channels include the e + jets and fi + jets subchannels. These events 
are divided into two disjoint subsamples: untagged (£ + jets) and (SLT) tagged (£+jets//i). 
The major background for both of these subsamples is from W + jets events. There is also 
a smaller component from QCD multijet events, in which one jet fakes a lepton signature. 

The cuts imposed for the untagged £ + jets channel are summarized in Table ||, the 
basic requirements being a high-p T central lepton, large $ T , four jets, and no 6-tag. For the 
first half of Run lb, the muon acceptance was restricted to \q\ < 1.0 because of radiation 
damage in the forward muon chambers. For Run la and the second half of Run lb (after 
the damage was repaired), the acceptance extended out to \rj\ < 1.7. The variable rj w shown 
in the table is the pseudorapidity of the jy-vector formed by combining the lepton and the 
$ T . The longitudinal component is found by constraining the lepton-neutrino pair to the 
W boson mass. A cut \r] w \ < 2.0 is imposed because D0 observes that VECBOS, which is 
used to model the kinematics of the W + jets background, does not agree well with data in 
the forward region.il The cut has very little effect on the top quark signal, however. The 
variable Ej, is defined as the scalar sum of the lepton momentum and the $ T . It is useful 
mainly for rejecting QCD multijets background; again, this cut does not have much effect 
on the tt signal. 

D0 also applies cuts in the (A, Ht) plane. The cut values are chosen through an ex- 
haustive search of possible cuts in order to optimize the expected precision of the cross 
section measurement. The search method used by D0, called the random grid search^ was 
a key development as it provides a highly efficient way to search n-dimensional parameter 
spaces. It differs from a conventional grid search in that the grid of cuts is determined by 
the expected distribution of the signal; hence the term "random" grid. The results of the 
search, along with the optimal cuts, are shown in Fig. The distribution of events in the 
(^4, Ht) plane is shown in Fig. |28|. Nineteen events survive all the cuts. 

The QCD multijets background is estimated from data by first measuring the probability 
for a jet to fake a lepton signature, and then directly applying these probabilities to the QCD 
multijets data. The dominant W + jets background is estimated using the technique of 
"Berends scaling." This is based on the observation that the number of iy + jets events falls 
off exponentially with the number of jets. Thus, after subtracting the contribution from 
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Fig. 27. Expected signal (for mj = 180 GeV/c 2 , e + jets channel) versus background yields in 
77 pb _1 for various cuts on A and Ht- The arrow identifies the optimal cut found by the random 
grid search method. 
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Fig. 28. Distributions the in (A,Ht) plane for data (top left), it MC (top right), multijet 
background (bottom left), and W + jets vecbos MC (bottom right). From Ref. 



QCD multijets background, D0 extrapolates from the number of £ + jets events observed 
at low multiplicities to estimate the number of W + jets events with at least four jets. 
The efficiency for these W + jets events to pass the (A, H T ) requirement is estimated using 
VECBOS. By this procedure, D0 estimates the background to be 8.7 ± 1.7 events. About 
15% of this background is due to QCD multijet events; the rest is from W + jets. 

D0 finds that the dominant systematic uncertainty of 15% is due to the differences in 
the (A, Ht) distributions between data and the corresponding simulations of W + 2 jet and 
iy + 3 jet events. The uncertainty due to possible departures from the Berends scaling law is 
found to be 10%. The jet energy scale uncertainty, however, is only 5%. With 19 candidate 
events and a background of 8.7 ± 1.7 events, the measured cross section for the t + jets 
channel is a ft (m t = 172.1 GeV/c 2 ) = 4.1 ± 2.1 pb. 

The cuts for the tagged channels i + jets//i are summarized in Table |]. Note that for 
these channels, only three jets are required, as opposed to four for the untagged channels. 
The (^4, Ht) cuts are also looser. This is to compensate for the loss of acceptance caused by 
requiring events to be tagged. To reject QCD multijets background, D0 makes the cuts:" 

• e+jets/> channel: $ T > 35 GeV if A0(# T ,//) < 25°. 

• n + jeta/fj, channel: A0($ T ,/i) < 170° and \A(j)($ T , /j.) - 90°|/90° < $ r /(45 GeV). 

In addition, in the /x + jets//x channel, D0 rejects Z — > background using a kinematic fit 
to the Z boson decay hypothesis. With these cuts, eleven candidates survive. 

The W + jets background is computed by starting with the assumption that the heavy 
flavor content of W + jets events is the same as that in multijet events. Therefore, from 
the tag probability per jet pj, as determined from a multijet sample (and parameterized in 
terms of Et and rj), D0 determines the tagged W + jets background by weighting a sample 
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Fig. 29. Jet multiplicity of 6-tagged lepton+jets events before A and Ht cuts. Square points 
represent data and the circles the expected background. From Ref. 



of background-dominated untagged events by the per-event tag probability J2jPj, where the 
sum runs over all jets. 

Figure |29| compares the observed numbers of events for different jet multiplicities to the 
background prediction. The agreement for the one and two jet samples is good, while for 
three or more jets there is a clear excess of events. The estimated background in the > 3 jet 
signal region is 2.4 ± 0.5 events, including a 10% systematic uncertainty due to the tagging 
procedure. The measured cross section for these channels is cr t j{m t = 172.1 GeV/c 2 ) = 
8.3 ±3.5 pb. 



5.3.3. All- Jets Channel 



The all-jets channel, as we noted earlier, suffers from a huge QCD background. It is 
a channel for which 6-tagging is particularly effective and D0, like CDF, takes advantage 
of this. However, since D0 cannot tag 6-jets nearly as efficiently as CDF, extensive use is 
made of multivariate methods, in order to exploit other, more subtle, differences between 
the signal and the background. The resulting analysisE3E2l is rather complex; therefore, we 
focus here on only the key ideas. 

The measurement of the ti production cross section in the all jets channel is based on 
an integrated luminosity of 110 pb -1 . The trigger used is sensitive to events with many jets 
and large Ht- The initial selection cuts are: 

• No isolated leptons. This ensures that there is no overlap with the other channels. 

• Ht > 115 GeV (using jets with R = 0.5, E T > 15 GeV, and \r]\ < 2.5). 

• At least six R = 0.3 cone jets and no more than eight R = 0.5 cone jets (with Et > 
10 GeV and \r)\ < 2.5). Both of these requirements improve the signal/background 
ratio: the lower bound rejects only 14% of the signal while removing 36% of the 
background, and the upper bound loses 5% of the signal but 13% of the background. 



59 



Table 10. List of variables used by D0 in its all-jets cross section measurement. The Qi are the 
eigenvalues (in increasing size) of the normalized momentum tensor = J2j PjaPjb/ J2j \Pj\ 2 of 
the jets in the event; a, b run over the x, y, z components of the jet momenta Pj. 



Variable 


Description 


Ht 


Total jet transverse energy 




Total tt center of mass energy 


Ht3 


Transverse energy of the non-leading jets 




Average jet count (see text for details) 


A 


Aplanarity (§Qi) 


S 


Sphericity (§(Qi + Q 2 )) 


C 


Centrahty (H T / £ iets ^ et ) 



Ej^ 1 /Ht Transverse energy fraction of leading jet 

\J E^ t5 E^ te Geometric mean of Et of the fifth and sixth jets 

(rj 2 ) Mean square rapidity of all jets 

p T Transverse momentum of the tag muon 

T Fisher discriminant: A measure of the difference 

in jet widths between signal and background 
Jv[ Mass likelihood: A measure of the degree to which 

jet invariant masses satisfy the constraints that there 

be two W bosons and that the two top quarks have 

equal mass. 



Additional cuts are applied to remove events with noise from the Main Ring. After these 
cuts, about 280,000 events remain. At this stage, the signal to background ratio is about 
1 : 1000. Requiring an SLT 6-tag increases the signal/background by about an order of 
magnitude, leaving 6000 events. 

To make further progress, D0 performs a multivariate analysis using thirteen variables, 
described briefly in Table [H]. The principal ones are Ht, Ht3 = Ht—E^—E^ 2 , the average 
jet count N^ ts , the aplanarity A, the centrality C = Ht/ Zljets E^ 1 , and the transverse 
momentum of the muon p T . A particularly powerful (and unusual) variable is the average 
jet count, denned by 

/■55 GeV / r 55 GeV 

N£ ts = / E T N(> E T ) dE T / / E T dE T , (17) 

J Jib GeV / J15 GeV 

where A^(> Et) is the number of jets with transverse energy greater than Et- This variable, 
inspired by the work of Tkachov,ffl is interesting in that it assigns a nonintegral "number of 
jets" to an event. 

These thirteen variables are combined using feed-forward neural networks into a single 
discriminant NN2, the value of which lies between zero and unity. This is done in three 
steps: 



- a network is trained with the first seven variables listed in Table 



using all 
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Fig. 30. (a) The mean muon pt versus jet Et, as observed by D0 in their all-jets cross section 
analysis. We see that the muon pt is independent of the jet Et- (b) The observed muon pt 
distribution for three jet Et ranges (normalized to the same number of events). From Ref. 89. 



events; this network is used to select better event samples for training the subsequent 
networks. 



iViVi — a network is trained with the first ten variables of Table [TO, having twenty 



hidden nodes and one output node, using events which satisfy the cut iViVo > 0.3. 

• iVA^ — a final network is trained, using events with iViVo > 0.3, with four input nodes, 
eight hidden nodes and one output node; the inputs are the output of NNi plus the 
remaining three variables. 

Since the data are predominantly background, D0 uses untagged events as the back- 
ground sample in the training, and tagged herwig it Monte Carlo events (with m t = 
180 GeV/c 2 ) for the signal. There is, however, a complication. Untagged events, by defini- 
tion, do not have muons; but the latter are nonetheless needed to train iVA^! The solution 
adopted is to generate Monte Carlo muons according to the measured muon pt spectrum 
and to add them to the untagged sample, thereby transforming the background sample into 
a tagged sample. One might worry that this procedure would fail because of correlations 
between the muon pt and the remaining variables. If correlations exist, the incoherent ad- 
dition of Monte Carlo muons to untagged events would not be expected to reproduce the 
characteristics of tagged events. 



However, the procedure works. As shown in Fig. |30], the muon pt is uncorrelated with 
the Et of the tagged jet, provided that the tagged jet is not corrected for the presence of 
the muon. This absence of correlation between muon pt and jet Et precludes the existence 
of correlation between the muon pt and the rest of the event. 
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Table 11. Summary of D0 cross section measurements, for mt = 172.1 GeV/c 2 . 



Channel 


e x BR(%) 


Data 


Background 


a tt (P b ) 


11 (with ev) 


0.91 ± 0.17 


9 


2.6 ± 0.6 


6.4 ± 3.3 


£ + jets (untagged) 


2.28 ±0.46 


19 


8.7 ± 1.7 


4.1 ±2.1 


i + jets //j, 


0.96 ±0.15 


11 


2.4 ±0.5 


8.3 ±3.5 


U and t + jets 


4.14 ±0.69 


39 


13.7 ±2.2 


5.6 ± 1.8 


All-jets (iViVo > 0.85, m t = 180 GeV/c 2 ) 


2.2 ±0.4 


41 


24 ±2.4 


6.5 ±3.0 


All-jets (from fit) 








7.1 ±3.2 


Total 








5.9 ± 1.7 



The background to the tagged data is estimated by weighting the untagged data (with 



the Monte Carlo muons) by the tag probability per event. Figure |3l] illustrates how well the 
weighting procedure works. The figure shows an absolute comparison of the tagged events 
to the predictions for the thirteen variables. For all variables, the agreement is seen to be 
good. 

To extract the cross section, D0 fits the observed iViV2 distribution for tagged events 
to the expectations for signal and background. The fit is performed for the 2207 tagged 
events that pass the cut NN2 > 0.02. The results of the fit (for m t = 180 GeV/c 2 ) are 
shown in Fig. [32[ This gives a cross section for m t = 180 GeV/c 2 of 6.3 ±2.5 (stat) pb. The 
major systematic uncertainties in this measurement are in the background tag probability 
(7%), the it tag probability (7%), the modeling of the muon px spectrum, and the jet 
energy scale (9%). After interpolating to the measured top quark mass, the final result is 
a ft (m t = 172.1 GeV/c 2 ) = 7.1 ± 2.8 (stat) ± 1.5 (syst) pb. 

As a cross check, D0 also estimates the cross section using the usual counting method. 
Here, a cut NN2 > 0.85 is imposed. This yields 41 observed events, with an estimated 
background of 24 ± 2.4 events. The signal efficiency at m t = 180 GeV/c 2 is 0.022 ± 0.0037, 
yielding a cross section of 6.5 ± 2.6 (stat) ±1.4 (syst) pb, in good agreement with that found 
by the fit. If the cut value is reoptimized to maximize the expected significance (using Monte 
Carlo), the result is NN 2 > 0.94, yielding 18 observed events and an estimated background 
of 6.9 ± 0.9 events. The chance of this background fluctuating to give at least the observed 
signal is 6 x 10~ 4 , or a 3.2cx deviation for a Gaussian. 

5.4. Summary of Cross Section Measurements 



The results of the CDF and D0 cross section analyses are summarized in Tables p| and [IT 
respectively, and graphically in Fig. [£| 

Excluding the all-jets channel, D0 observes a total of 39 events with a background of 
13.7 ±2.2 events. For a top quark mass of 172.1 GeV/c 2 this gives a measured cross section 
of 5.6 ± 1.8 pb. If the measurement from the all-jets channel is included, the result changes 
to 5.9 ± 1.7 pb. The measurement from CDF, which includes all channels except the r- 
dilepton channels, is a t j = 7.6ll\[l pb for a top quark mass of 175 GeV/c 2 . Both results are 
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Fig. 31. A comparison of D0's absolute predictions (histograms) for the distributions of the 
thirteen variables used in the all-jets cross section measurement: (a) Ht (GeV), (b) y/§ (GeV), (c) 



4 etl /#T, (d) Ht3 (GeV), (e) JV*,, (f) yJl^B^ (GeV), (g) A, (h) S, (i) C, (j) (r? 2 ), (k) p£, (1) 
J 7 , and (m) A4. The points represent the observed data in a sample of 3853 tagged events. From 
Ref. M 
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Fig. 32. The distribution of the final neural network output (NN2) in the D0 all-jets cross 
section analysis showing the results of a fit of the observed tagged distribution to the predictions 
for signal and background, with nit = 180 GeV/c 2 . Only statistical errors are shown, and events 
with NN2 < 0.02 are not plotted. The fit has a x 2 P er degree of freedom of 16.9/17. From Ref. 88. 
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presented graphically in Fig. [34], which compares the measurements to several predictions. 
Although the cross sections obtained by CDF are typically higher than those from D0, both 
sets of results are consistent with each other and with predicted values. The results are also 
consistent with the expected branching ratios. 



6. MEASUREMENT OF THE TOP QUARK MASS 

Early direct measurements of the top quark mass were published in the papers announc- 
ing the observation of the top quark in 1995. @'iH Since that discovery, the two collaborations 
have mounted prodigious efforts to measure the top quark mass as precisely as possible. It 
is to these results that we now turn. After sketching the analysis strategies used by the two 
collaborations, we describe in greater detail interesting aspects of the analyses of CDF and 
D0. We end with a summary of the results. 



6.1. General Strategy 

Both collaborations use two different strategies to measure the mass of the top quark. 
The first involves kinematically fitting the events to a tt hypothesis and then extracting 
from the fits an estimate of the mass of the top quark. The second strategy is based on 
constructing mass- dependent variables that do not require a kinematic fit. 

In the first strategy, one tries to identify the observed objects (jets, leptons, and missing 
transverse energy) with those at the partonic level. To the degree that such an identification 
is possible, one can perform a kinematic fit of the observed event to the tt decay hypothesis. 
When the mass of the top quark is extracted from such a fit, it is commonly called the fitted 
mass. Because the tt decay hypothesis does not always hold true, this is not the same as 
the true top quark mass; nevertheless, it is a well-defined observable, the mean of which 
depends on the mass of the top quark. This defines one class of mass- dependent observables 
or variables. 

The principal advantage of a kinematic fit is that, generally, the resulting variables have 
the greatest sensitivity to the mass of the top quark. But such fits are not always possible; 
this is so, for example, if there are fewer observed jets than partons. 

The alternative strategy eschews the complexities of a kinematic fit, and instead uses 
one or more kinematic quantities that are known to be mass-dependent. An example of such 
a variable is = YIE^ , the scalar sum of the jet transverse energies. The disadvantage 
of this approach is that variables such as Hj- tend to have much broader distributions than 
those of the fitted mass. On the other hand, this approach is simple and direct (compared 
with kinematic fitting), and it can be used for all events. Moreover, it may ultimately 
be the more satisfactory strategy for the following reason. As the accuracy of the mass 
measurements improve, we shall reach a point when the quantum nature of partonic decay 
can no longer be neglected and our classical one-to-one identification of a jet with a parton 
will break down.0 When this point is reached, it may be necessary to use observables that 
do not require forcing a one-to-one map back to the partonic level. Ultimately, the limitation 
will be the sensitivity of the mass to the QCD modeling of the final state. 
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Fig. 33. The tt production cross section measured in the channels studied by CDF and D0. 
Also shown is the range of predictions from various theoretical calculations. 
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Fig. 34. The tt production cross sections as measured by CDF and D0. Also shown are various 
theoretical calculations. 

Each strategy furnishes a mass-dependent variable x per event, the distribution of which, 
f(x), depends upon the mass of the top quark m t if the event is a top quark event. The func- 
tions f 8 {x\m t ) and /&(#) pertaining to the signal and background, respectively, are calculated 
using detailed Monte Carlo simulations of the signal and background reactions. The func- 
tions f s (x\mt) and fb{x) are sometimes referred to as templates. Given a set of templates and 
a set of measurements (x±, ■ ■ • , Xn), there are a variety of standard methods available, most 
commonly maximum likelihood^ to extract a single estimate of the parameter of interest, 
here the mass m t . 

Much of the effort in these analyses is devoted to constructing reliable templates. This is 
a delicate matter because of the necessary reliance on complicated Monte Carlo simulations 
to compute them. The key issue, of course, is to establish that the templates, so obtained, 
agree with the data. For the background, this is less of an issue because the data, being 
largely background, afford a direct check of the correctness of those templates in regions of 
parameter space that are relatively free of signal. However, for the signal, the verification of 
the templates is, necessarily, more problematic because the signal events are generally less 
numerous. It is therefore notable that in spite of these difficulties, both collaborations have 
been able to produce very precise and consistent results. 

At a fundamental level, there is really no difference between the two strategies we have 
outlined: the result of the kinematic fit is a mass- dependent kinematic quantity, albeit a 
particularly complicated one. Nevertheless, the practical details of carrying out the analyses 
are sufficiently different that the distinction is worth making. 

In the following sections we describe how these strategies have been realized by each 
collaboration. We begin with the analyses of CDF, followed by those of D0. 
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6.2. CDF Analyses 

The CDF collaboration has measured the top quark mass in all of the dilepton, lep- 
ton+jets, and all jets channels, using their full Run 1 data sample. 



6.2.1. Dilepton Channels 

In this analysis,!^ CDF considers the channels e/i, /i/i, and ee, in which both W bosons 
decay leptonically. The mass analysis starts with the same data sample and cuts as used 
in the dilepton cross section analysis (see Sec. |57^ and Table ||). To improve the sig- 



nal/background ratio beyond the level used for the cross section analysis, CDF imposes 
the cut Hf > 170 GeV, where Hf = p 1 ^ + p^ + E^ 1 + E$ t2 + $ T is the sum of the 
transverse momenta of the two leptons, the transverse energies of the two highest-.&r jets, 
and the missing transverse energy. The power of such variables to discriminate signal from 
background was recognized early onjl and was first applied to good effect by D0.E1EI With 
this final cut, the sample is reduced to eight events with an estimated background of 1.3±0.3 
events. The CDF collaboration then estimates the top quark mass using two different meth- 
ods. 

Both methods are a realization of the second strategy, being based on relatively simple 
mass-dependent kinematic variables. For the first method, CDF uses the observation that 
the energies of the 6-quarks are sensitive to the mass of the top quark, m t . Indeed, it is found 
that the mean energy of the two highest jets depends linearly on m t , with a slope of about 
0.5. From Monte Carlo simulations, CDF obtains a set of templates f(E\m t ) describing the 
observed jet energy distribution, one for each top quark mass considered in the range 100 
to 240 GeV/c 2 . Templates are also constructed for the backgrounds. CDF then performs 
a maximum likelihood fit of the observed distribution of jet energies to a mixture of ti and 
background events, with the background count constrained to 1.3 ± 0.3 events. For each 
assumed top quark mass m, — In (C(m)), where C(m) is the value of the likelihood at that 
mass, is calculated. By fitting a third order polynomial to the negative logarithms, and using 
the position of its minimum as an estimate of the mass, CDF obtains the mass estimate 
159 ±23 (stat) ±11 (syst) GeV/c 2 . Following standard practice, CDF defines the statistical 
error by applying the "0.5-rule": the error is the mass difference, relative to the value at 
the minimum (159 GeV/c 2 ), that increases — ln(£) by 0.5. Figure 35a shows the jet energy 
distribution and the resulting log-likelihood function. 

The second method exploits the relationship between m t and the invariant mass 
of the lepton and the associated 6-quark. Since the top quark decays to a 6-quark and 
a W boson, we can write mf = m 2 ± + 2E b M w , where E b is the energy of the b- 
quark in the rest frame of the W boson. The 6-quark's energy E b , which is constant in the 
W boson rest frame, can be calculated from the mean invariant mass (M#,) and opening 
angle (cos 6^) between the 6-quark and the lepton: (Mf b ) = E b M w {\ — (cos 6^)). To 
lowest order in the Standard Model, (cos 6^) = M^/(m 2 ± 2M^), which leads to m 2 = 

(M| 6 ) + \Jm^ + 4M^ (M| b ) ± (M| 6 ) 2 . Given this relationship, a measurement of (Mjj,) is 
tantamount to a direct measurement of m t . (For an interesting variation on this theme, see 
Ref. H) 
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Fig. 35. (a) The energy distribution of the two highest Et jets for the CDF dilepton events 
(histogram, two entries per event). The lightly shaded curve is the prediction for it signal 
(m t = 160 GeV/c 2 ) plus background, and the heavily shaded curve is the prediction for back- 



ground alone. The inset shows the — ln(£) fit as a function of rat- (b) The distribution of Mj 



(histogram, two entries per event). The shaded curves are as in (a). From Ref. 81 
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Unfortunately, there is a difficulty: a priori, we do not know which jet should be paired 
with which lepton; whether (£ + ji,£~j2) is the correct (£ + b,£~b) assignment or (£ + j2,£~ji). 
In addition, since is defined by a quadratic equation, there are two solutions to be 
considered for each jet configuration. CDF resolves the first of these difficulties by selecting 
the jet configuration with the smaller sum of invariant masses. This selection is correct 
55% to 75% of the time, depending upon the top quark mass. The values of Ma, selected 
in this manner are denoted M^ min . The second difficulty is dealt with by retaining both 
solutions, and plotting each event twice. Figure [35| b shows the resulting M^ min distribution 
for the eight events in the sample, along with the signal and background predictions. Of 
course, owing to selection biases and jet energy mismeasurements, (M^^J) is not the same 
as (Mfb); however, from Monte Carlo studies, CDF finds that the two quantities are linearly 
related. The result from this method is m t = 163 ± 20 (stat) ± 9 (syst) GeV/c 2 . Combining 
this with the previous result, taking into account their mutual correlation, gives m t = 
161 ± 17 (stat) ± 10 (syst) GeV/c 2 . 

For both methods, the major sources of systematic uncertainties are the knowledge of 
the jet energy scale and the shape of the background distributions. 

At a recent conference in Vancouver (ICHEP98), CDF presented an updated dilep- 
ton channel mass measurement,^ using a method similar to the D0 z/WT method (see 
Sec. |6.3.1| ). The result of this analysis was m t = 167.4 ± 10.3 (stat) ± 4.8 (syst) GeV/c 2 . 



(The details of this analysis were not available in time for inclusion in this review.) 



6.2.2. Lepton+Jets Channels 

These channels provide CDF with its most precise measurement of the top quark massif! 
This is due to both the relatively large branching ratios for these channels (~ 30%) and the 
large signal/background ratio that can be achieved through 6-tagging (see Sec. |3.3.6|) . The 
event selection used is similar to that used for the cross section analysis (see Table f|), 
except that four jets are required, and events without a 6-tag are accepted. But in order to 
increase the efficiency for the four jet requirement, the cuts on the fourth jet are loosened to 
E T > 8 GeV and \r}\ < 2.4, provided that one of the four leading jets is 6-tagged with either 
the SLT or the SVX algorithm. (For untagged events, all four jets must satisfy Et > 15 GeV 
and \rj\ < 2.0.) This selection yields a sample of 83 events. 

This analysis follows the strategy of kinematic fitting. The observed event is fit to the 
hypothesis pp — > it + X — > (£ub)(bqq) + X, with the four leading jets assumed to map onto 
the four quarks. The measured variables are the three-momenta of the lepton and the four 
leading jets, plus the vector sum of all remaining transverse energy (excluding the lepton 
and the leading jets). This latter vector gives the transverse momentum of the recoiling 
system X. The momentum of the neutrino is not measured, yielding three unknowns. The 
neutrino and lepton are assumed to be massless, 6-quarks are assigned a mass of 5 GeV/c 2 , 
and light quarks are assigned a mass of 0.5 GeV/c 2 . Five constraints can then be imposed: 

• The transverse momentum components of the it + X system are zero. 

• The invariant mass of the lepton and neutrino equals M w . 
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• The invariant mass of the q and the q equals M\y- 

• The mass of the t equals that of the t. 

Since there are two more constraints than unknowns, this is a 2C fit. It is solved by standard 
X 2 -minimization techniques, resulting, for each event, in a fitted mass mat and a x 2 value 
that is taken as a measure of how well the event satisfies the it hypothesis. 

The energy of each jet is adjusted to match as closely as possible the energy of the 
parton from which the jet is presumed to have originated,^ starting with corrections of the 



sort outlined in Sec. |3.3.1| . In addition, CDF applies energy corrections derived by directly 
comparing, in herwig Monte Carlo samples, the energies of partons with the reconstructed 
jets. These corrections depend on parton flavor; hadronically decaying 6-quarks, leptonically 
decaying 6-quarks, and light quarks are all corrected differently. Clearly, these corrections 
depend critically on the ability of herwig to model the fragmentation properties of jets. It is 
therefore pertinent to ask how well these properties are modeled. The answer is, surprisingly 
well, given the phenomenological nature of the fragmentation model. Detailed studies of the 
energy flow within jets (see, e.g., Ref. |96|) confirm the ability of herwig to model jets 
accurately. 

In the lepton+jets channels, there is considerably more ambiguity in assigning the jets 
to the partons than is the case for the dilepton channels. There are twelve possible config- 
urations, or assignments of the four leading jets to the four partons. Moreover, one must 
choose an initial value for the longitudinal component of the neutrino momentum. Using 
one of the constraints, this quantity can be found to within a two-fold ambiguity. To guard 
against finding a local minimum, the fit is tried with both starting points, giving a total of 
up to twenty-four kinematic fits per event. However, jets that are tagged are only assigned 
to 6-quarks, thus reducing the number of configurations. The configuration with the small- 
est x 2 is chosen, but events whose smallest x 2 value exceeds 10 are rejected. This x 2 cu t 
reduces the sample from 83 to 76 events. 

The x 2 variable, unfortunately, is not very effective at selecting the correct configuration 
(see the discussion in Sec. |6.3.2|) . Consequently, although the standard deviation of the m fit 
distribution would be only ~ 13 GeV/c 2 (for a top quark mass of 175 GeV/c 2 ) were the 
selection procedure perfect, CDF finds that, in practice, the standard deviation of the m fit 
distribution is approximately double that value. 

The CDF collaboration uses a maximum-likelihood method to extract an estimate of the 
top quark mass from the sample of 76 events. The likelihood function has the form 

C = C shape (D\m t ,x b ,a) x £backgr(z&k°) x C paTauia (a\a ) , (18) 



"Note, however, that this procedure is not strictly necessary provided that the Monte Carlo 
simulation offers a faithful rendition of the characteristics of the data, the most critical being that 
the energy scale of the simulation agrees with that of the data. Once this has been achieved, the 
only requirement is that the same kinematic fitting algorithm be applied both to real and Monte 
Carlo events. The rescaling may, however, improve the resolution slightly because of the imposition 
of the Myy constraints. 
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where 

N 



C Bhape (D\m t ,x b ,a) = JJ[(1 - x b )f s {mi\m tl a) + x b f b (mi\a)] (19) 



i=i 



is the joint probability density assigned to the sample of N fitted masses D = (mi, . . . , mjv); 
f s (mi\m t , a) is the probability density to observe the fitted mass rrii given a top quark mass 
m t , assuming that the ith event is signal and f b {rrii\a) is the probability density to observe 
rrii assuming the event is background; x b is the background fraction; and a denotes the set 
of parameters that determine the shape of the functions f s (rrii\m t , a) and f b (jrii\a). These 
functions (templates) are derived from Monte Carlo calculations (based on herwig for top 
events and vecbos for the W+jets background). In previous work,0 CDF used a discrete 
set of templates. However, the collaboration now finds that smaller systematic uncertainties 
result from deriving from the discrete set a smooth parameterization of / s (mj|mt, a), both 
in rrii and in rrit- (Amusingly, D0, which earlier used parameterized templates, now finds it 
better to use a discrete set!) 

These functions contain parameters a, other than m t , that define the shape of the func- 
tions. The parameters have estimates a and uncertainties that are encoded in the likelihood 
function £ para m- This likelihood function constrains the shape of the templates in £ s h ap e to 
vary within their uncertainties. The shape parameter uncertainties reflect simply the finite 
statistics of the Monte Carlo event samples that were used to generate the discrete set of 
templates. 

The background fraction x b is constrained within its uncertainty to the independently 
measured background fraction x® by the background likelihood function £backgrOE&|#&) • One 
can think of £back g r and £ param as prior probabilities with respect to the likelihood function 
£shape- They encode the information (a°,x°) that is known about the background fraction 
and the shape of the templates. 

In statistical analysis, it is sometimes advantageous to use stratified sampling;!! that is, 
to divide a sample into two or more subsamples. This is especially true when the character- 
istics vary widely from one subsample to another. In order to make optimal use of all the 
available information, the event sample is divided into disjoint subsamples. Each subsample 
is assigned its own likelihood function in accordance with Eq. fll8|), the product of which is 
then maximized to obtain the final estimate of the top quark mass and its uncertainty. From 
Monte Carlo studies, CDF concludes that four subsamples give the optimum partition: 

• events with a single SVX tag, 

• events with two SVX tags, 

• events with an SLT tag but no SVX tag, and 

• events with no tag but with all four leading jets satisfying Et > 15 GeV and \rj\ < 2. 

The expected background fraction x® in each subsample is derived from the W+ > 3 jets 
background estimate used in the it cross section measurement (see Sec. |5.2| ). To estimate the 
background in the subsamples, one needs the efficiencies of the additional cuts, including the 
fourth jet requirement and the x 2 cut. These efficiencies are determined from Monte Carlo 
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Table 12. Subsamples used in the CDF lepton+jets top quark mass measurements. N Q \> S is the 
the number of observed events, x\, is the expected background fraction, and mt is the measured 



top quark mass. Uncertainties in mt are statistical only. From Ref. |95 







x b 


Measured mt 


Subsample 


N oha 


(%) 


(GeV/c 2 ) 


SVX double tag 


5 


5±3 


170.1 ±9.3 


SVX single tag 


15 


13 ±5 


178.0 ± 7.9 


SLT tag (no SVX) 


14 


40 ±9 


142.0ifl 


No tag 


42 


56 ± 15 


181.0 ±9.0 



studies. Using the tagging efficiencies from the cross section analysis and the background 
rates for the subsample, CDF writes the total mean number of events n in each subsample 
as a function of the mean number of signal and background events S and B, respectively, in 
the combined sample; that is, n = e s S + ebB, where e s and e& are known fractions, derived 
from the Monte Carlo studies. The parameters S and B are then estimated by maximizing 
a multinomial likelihood function, based on the measured event count Cj in subsample i and 



its mean count n,. This leads to the background fractions = ebB /m, listed in Table 12 



It is found that about 67% of the background is due to W + jets, 20% is due to multijet 
events (in which a jet has been misidentified as a lepton) and bb events with the 6-hadron 
decaying semileptonically, and the remaining 13% comes from [Z — > ££)±jets events, diboson 
events (WW, WZ, and ZZ), and single top. 

As noted above, a key issue is the construction of reliable templates. For the signal 
templates, one can judge their reliability only after the mass distributions have been fitted. 
If the fit is good, this provides ex post facto evidence of their soundness. The reliability 
of the background templates is generally addressed by comparing Monte Carlo simulations 
with data in regions that are largely devoid of signal. If the templates agree well in those 
regions, then one has some confidence that they are probably adequate in the signal-enhanced 
regions. These cross checks are not rigorous: the fact that the templates are modeled well 
in the background-dominated regions does not guarantee that they are also modeled well in 
the regions of interest. This observation applies, of course, not only to mass templates but 
to any distribution based on Monte Carlo simulation. However, the quoted uncertainties 
account for possible mismodeling of the templates. Both CDF and D0 confirm the reliability 
of the templates in this manner. 

Each subsample provides an independent mass measurement, as listed in Table [12. 
The fitted mass distributions and corresponding log-likelihood curves for each are shown 
in Fig. |36|. The combined mass distribution, along with the total log-likelihood, is shown in 
Fig. [IT]. From this likelihood curve, CDF measures a top quark mass of 175.9 ±4.8 GeV/c 2 , 
where the uncertainty is defined by the usual "0.5-rule" (described above). The dominant 
systematic uncertainty, of 4.4 GeV/c 2 , is from the jet energy scale; see Table [TBI. A very 
important goal for the next round of experiments is to substantially reduce this uncertainty. 
This will be no easy task; an enormous effort has already been expended to achieve the 
present small uncertainty. 
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Fig. 36. The CDF top quark mass distributions for the four lepton+jets subsamples. The 
points are the data, the heavily shaded area is the top signal + background combination resulting 
from the fit, and the lightly shaded area is background alone. The insets show the negative log 



likelihood vs. m< for each subset. From Ref. 95 
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Fig. 37. The CDF top quark mass distribution for all four subsamples combined. The points 
are the data, the heavily shaded area is the top signal + background combination resulting from 
the fit, and the lightly shaded area is background alone. The inset shows the negative log likelihood 
vs. rat. From Ref. 95. 



75 



Table 13. Systematic uncertainties in CDF's rat measurement in the lepton + jets channels. 



From Ref. |j. 

Source Value (GeV/c 2 ) 

Jet energy measurement 4.4 

Initial and final state radiation 1.8 

Shape of background spectrum 1.3 

b tag bias 0.4 

Parton distribution functions 0.3 

Total 4~9 



6.2.3. All Jets Channel 

CDF has also presented a top quark mass measurement in the challenging all-jets 
channel,il in which both W bosons decay hadronically. The event sample is based on the 
kinematic cuts of technique I of the cross section measurement (see Sec. |5.2.3| ), except that 
six jets are required and the cut on Ht = J2 is relaxed from 300 GeV to Ht > 200 GeV. 
The latter change is to reduce the bias in the mass distributions due to threshold effects. 
Before requiring a 6-tag, 1121 events survive these cuts. After tagging, 136 events remain. 

CDF obtains a fitted mass for each event using a kinematic fit to the hypothesis it — > 
(bW + )(bW~) — ► {bqq)(bqq). One can parameterize these events so that there are thirteen 
unknowns (the top quark mass and the three-momenta of both top quarks and W bosons) 
and sixteen constraints (four-momentum conservation at both top quark and W boson decay 
vertices). Thus, a 3C fit is possible. As in the kinematic fit in the lepton + jets channel, 
CDF tries all possible assignments of jets to the quarks which are consistent with assigning 
the tagged jet as a 6-jet. The fit which gives the smallest \ 2 is retained. Figure [38| shows 
the resulting three-jet fitted mass distribution. 

A major challenge for this analysis is the background estimation. This is because not 
only is the background large, but its fitted mass distribution is not that different in shape 



from the signal (see Fig. [38]). Therefore, an error in the background normalization could 
result in a bias in the mass measurement. CDF estimates the background using a technique 
similar to that used for the lepton+jets and all-jets cross section measurements: each jet in 
the sample of 1121 untagged events is weighted by the probability for a jet in multijet events 
to be tagged. This tag probability^ is parameterized in terms of the jet E T , rj, and track 
multiplicity and the event aplanarity. This technique yields an estimated background of 



The appropriate parameterization of the tag probability has been a source of seemingly endless 
discussion within CDF and D0. Several different parameterizations have been tried, whose forms 
are largely the outcome of a judicious use of intuition and experimentation. The chief difficulty 
is to understand whether the choice of parameters, in terms of which the tag probability is pa- 
rameterized, is complete in the sense that they provide a full account of the variation of the tag 
probability from event to event. Yet, in spite of this and other difficulties, the tag probability 
functions work remarkably well in practice. 
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108 ± 9 events. The shape of the fitted mass distribution is taken directly from the untagged 
sample. 

The final mass estimate is extracted from the fitted mass distribution using a maximum 
likelihood fit based on discrete templates (similar to what was used in the discovery paperB). 
The result is 186 ± 10 GeV/c 2 (see Fig. |38| ). In the published reference,^ CDF quotes a 
systematic uncertainty due to the modeling of gluon radiation and fragmentation effects of 
±8.6 GeV/c 2 . However, this uncertainty was recently revised downwardsS to ±3.0 GeV/c 2 . 
The uncertainty due to the jet energy scale is ±5.4 GeV/c 2 . Other systematic uncertain- 
ties are due to the fitting method (±5.2 GeV/c 2 ) and the background estimation method 
(±1.7 GeV/c 2 ). When these are combined in quadrature, one obtains a total systematic 
uncertainty of ±8.2 GeV/c 2 . 



6.3. D0 Analyses 



The D0 collaboration has measured the top quark mass in both the dilepton and lep- 
ton±jets channels. 



6.3.1. Dilepton Channel 

A dilepton it event is kinematically underconstrained due to the two unobserved neutri- 
nos. This would seem to preclude the use of kinematic fits as a means to estimate the top 
quark mass. However, Kondoi and Dalitz and GoldsteirE have shown that it is possible to 
estimate the mass from an ensemble of kinematic fits by making use of additional information 
concerning the likelihood of given final states. These suggestions were pursued vigorously 
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by D0 and led to a successful measurement of the top quark mass in these channels. E2 

The D0 collaboration has five dilepton candidate events (3 e/i, 1 ee, and 1 fifi) after 
applying the selection criteria used in the cross section measurement. An additional ee 
event is selected by relaxing the track requirement on one of the electrons in the presence 
of a semileptonic (soft muon) b-t&g. These six events are used to measure the top quark 



mass. 
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The expected backgrounds are 0.21 ± 0.16, 0.47 ± 0.09, and 0.73 ± 0.25 events 
in the e/i, ee, and /x/i channels, respectively. 

The invariant mass constraints m(£ii>i) = m^v-i) — an d m{l\V\b\) = m(£ 2 ^2^2) 
depend on the unknown neutrino four- vectors, v\ and z/ 2 . But, given a value for the top 
quark mass m t , the system can be solved to within a four-fold ambiguity. The crucial 
observation is that not all hypotheses about the value of the top quark mass are equally 
probable: some final-state kinematic configurations are more likely than others, and so some 
top quark masses are preferred over others. Therefore, if one could assign such probabilities 
given the measured momenta, it would be possible to estimate the top quark mass. For an 
assumed value of m if D0 assigns a probability p(D\m t ) to the measured momenta D and 
uses Bayes' theorem. 
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p(D\m t )p(m t ) 
/ °° p(D\m' t )p(m' t ) dm' t 



PKP) = ^=^-7 , (20) 
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Fig. 38. The CDF fitted mass distribution for the all jets sample. The points are the data, 
the shaded histogram is the background, and the open histogram is the sum of it signal with 
rat = 175 GeV/c 2 and background. The inset shows the log-likelihood curve. From Ref. p5|. 
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to compute the probability p(m t \D) of m t given the fourteen measured parameters compris- 
ing the data D. The expression p(m t ) is the prior probability for m t . 

For each value of mt, there are up to four possible solutions per jet-to-parton assignment, 
or configuration. For each solution and configuration, a likelihood L is computed. These 
likelihoods are summed over all solutions and configurations to arrive at the overall likelihood 
p(D\m t ). D0 chooses the prior probability for m t to be flat, so that Eq. (|20|) reduces to 
p(m t \D) oc p(D\m t ). 

In principle, the likelihood can be computed analytically, using 



Here, v is the set of eighteen parameters needed to completely specify the kinematics of 
the event, M. is the matrix element for the process qq — > tt + X — > £ + ub£~ub, and f(x) 
and f(x) are the parton densities for quarks and antiquarks of momentum fractions x and 
x in the proton and antiproton, respectively. The detector resolution function r(D\v) gives 
the probability to observe the data D given the true kinematic configuration v. The four- 
dimensional delta function 5 4 enforces the kinematic constraints: 



S 4 = 5(m(i + v) - M w )6(m(£~9) - M w )5(m(i + vb) - m t )5{m{£-pb) - m t ). (22) 



Unfortunately, it is not feasible to evaluate this expression for the large number of Monte 
Carlo events which are used to construct the templates. Therefore, D0 computes only much 
simplified approximations to this full likelihood. Two different methods are employed for 
doing so, as discussed below. 

The first method, called matrix element weighting (AiWT), is a modified version of the 
procedure suggested by Dalitz and Goldstein.ta Here, the sum of the transverse momenta 
of the two neutrinos is required to be equal to the measured $ T . The system is then solved 
for the neutrino momenta, and hence the top and antitop four-vectors. The likelihood for a 
given solution and configuration j is taken to be 



where f(x) and fix) are again the parton distribution functions (D0 uses CTEQ3ME3) and 
P{Ef M \rrit) is the probability density for the lepton energy in the rest frame of the top 
quark, Ef M . 

In the second method, called the neutrino weighting (z/WT) method, for each configura- 
tion j, the expected phase space of neutrino pseudorapidity in tt events (at a given m t ) is 
divided into elements with equal phase space weight. For each pair of neutrino rj values, a 
solution is sought. If a solution is found, the likelihood L(D\m t ,j) is assigned based on the 
degree to which the sum of the neutrino transverse momenta agrees with the Ip T in the event. 
For this calculation, only the smearing due to the underlying event is taken into account. 
For both methods, one sums L over all solutions and configurations to obtain p(D\m t ), and 
hence p(m t \D). 

So far, we have not accounted for the detector resolution function, r(D\d), where d are 
the true values of the momenta. (The information in d is a subset of that in v.) To include 
the effects of r, D0 uses the true, rather than the measured, momenta to solve the system. 




(21) 



L{D\m t ,j) oc f{x)f(x)P(Ef 1 M \m t )P{Eg M \m t ) 



(23) 
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But since the true momenta are unknown, it is necessary to consider all possible values of 
the true momenta d, constrained (probabilistically) by the measured momenta D, which of 
course are fixed. (In Ref. |100| , it is stated that the measured momenta are "smeared," but 
this we regard as physicists' argot for the statement we have just made.) The measured 
momenta D and the true momenta d are related via the experimental resolution function 
r(D\d), which modifies the relationship between p(D\m t ) and L(d\m t ,j) (now written in 
terms of the true momenta d) thus: 

p(D\m t ) cc I r(D\d)J2L(d\m t ,j), (24) 

j 

where the integration is over all possible momenta. D0 models its jet resolution function by 
a double Gaussian, which accounts for the inherent energy resolution of the D0 calorimeter 
plus the effects of large angle gluon radiation. The ]£ T is computed from the true momenta 
plus the contribution from the underlying event. The latter is included by fluctuating each 
component of $ T with a Gaussian, with a standard deviation of 4 GeV. 

Because of initial or final state radiation (ISR or FSR), an event may have more than 
two jets. In that case, D0 considers all possible interpretations of the jets, merging jets 
classified as FSR with the appropriate 6-jet. Each interpretation is assigned a weight, which 
is the product of the weights assigned to each jet. Therefore, for events with three or more 
jets, the sum in Eq. ( p4|) becomes a weighted sum. The weight assigned to an ISR jet is 
exp(— Et sin 9/(25 GeV)), where 6 is the polar angle of the ISR jet. For an FSR jet coming 
from a 6-quark, which together with the 6-jet form an invariant mass m, the weight assigned 
is exp(-m/(20 GeV/c 2 )). 

The normalized distributions of p(D\m t ) for the six candidate events are shown in Fig. RP. 



If the background in this sample of six events were negligible and no approximations were 
made in computing p(D\m t ), one could simply form the product of the six distributions and 
use the position of its peak or its mean as an estimate of the top quark mass. However, 
the background cannot be neglected, and the effects of the approximations must be taken 
into account. Therefore, D0 performs a maximum likelihood fit to a sum of signal and 
background templates. The signal templates are constructed from Monte Carlo simulations, 
the background templates from a combination of Monte Carlo and data. 

To make use of the information contained in the shape of the distributions p(m t \D), 
the latter are normalized and divided into five bins of width 40 GeV/c 2 (see Fig. |4"0| ). For 
each event i, the integrated content of each of the first four bins form the components of a 
four- dimensional vector Wi. Only four bins are used because the content of the fifth bin is 
determined by the normalization condition, and is therefore not independent. The likelihood 
function is then 

rt-\ \ n(-\ \ T3( at\ i ,-^nsfs(wi\m t ) +n b f b (wi) 

C(w\m t , n s , n b ) = G(n b \n b ) P(N\n s + n b ) I I ■ , (25) 

\ n s + n b 

where n s and n b are the mean signal and background counts, G(n b \n b ) is a Gaussian con- 
straint on n b , P(N\n s + n b ) is a Poisson constraint on (n s + n b ) to the observed sample size 
iV = 6, and f s and f b are four-dimensional probability density templates for the signal and 
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Fig. 39. Distributions of p(D\mt) for the six D0 dilepton candidates, using the AiWT (dashed) 
and i/WT (solid) methods. From Ref. |100| . 
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Fig. 40. The likelihood curve for a typical Monte Carlo ti event, normalized to unity. The 



vertical lines show the five bins into which the likelihood is divided. From Ref. 101 



Table 14. Summary of D0's dilepton mt measurements for full and partial data sets. Uncer- 



tainties are statistical only. From Ref. 100 



Channels Fit 


MWT (GeV/c 2 ) 


uWT (GeV/c 2 ) 


efj, + ee + 


168.2 ± 12.4 


170.0 ± 14.8 


e\x + ee 


168.0 ± 12.7 


173.3 ± 14.0 




173.1 ± 13.3 


170.1 ± 14.5 



The 



background, respectively. These templates are approximated using a multivariate probabil- 
ity density methodic The maximum likelihood estimate of m t and its error are determined 
by a quadratic fit to the — ln£ curve, using nine points about the minimum. 

The sums of the weights for the six candidate events in the five mass bins are com- 
pared with the signal and background expectations in Fig. (4~T , 
— In £ distributions from which the estimates m t 
m t = 170.0 ± 14.8 GeV/c 2 (z/WT) are extracted. Table [L4| shows the results of fits to sub- 
samples of the data. As shown in Table |15, the total systematic uncertainty is estimated 
to be 3.6 GeV/c 2 , dominated by the jet energy scale uncertainty. D0 combines the results 
from the AiWT and z/WT analyses, taking into account the 77% correlation between them, 
to obtain m t = 168.4 ± 12.3 (stat) ± 3.6 (syst) GeV/c 2 . 



insets show the 
168.2 ± 12.4 GeV/c 2 (MWT) and 



6.3.2. Lepton+Jets Channels 

The distributions of fitted mass (m^ t ) for signal and background can overlap significantly, 
depending on the mass of the top quark. Therefore, in order to extract the top quark 
mass reliably from an mg t distribution, it is necessary either to suppress the background 
sufficiently or to make optimal use of available information. Requiring the presence of fe-jets, 
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Fig. 41. The sum of the normalized candidate likelihoods for the .MWT (left) and z/WT (right) 
analyses, compared to the results of the maximum likelihood fit. The solid histogram is the best-fit 
it signal plus background, and the dotted histogram is background alone. The insets show the 



ln£ as a function of mt- From Ref. 100 



Table 15. Systematic errors for D0's measurement of mt in the dilepton channels. From 
Ref. BM 



Source Error (GeV/c 2 

Jet energy scale 2.4 

Signal model 1.8 

Multiple interactions 1.3 

Background model 1.1 

Likelihood fit 1.1 

Total 3.6 
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as we have seen, is a simple and effective means of enhancing the signal to background ratio. 
However, D0 uses only soft-lepton tagging, for which the tagging efficiency is only « 20% 
in lepton+ > 4 jets events, compared to ~ 53% at CDF, which has the ability to tag 6-jets 
with the SVX. 

It is therefore noteworthy that in spite of D0's relatively poorer 6-taggineJt has managed 
to measure the top quark mass with a precision approaching that of CDF.lM2a D0 achieves 
this by using multivariate techniquesBi!'^ for separating signal and background while 
minimizing the correlation of the selection with the top quark mass. Two multivariate 
methods are used to compute, approximately, a signal probability p(top|D) for each event, 
given data D. A likelihood fit of the data to a discrete set of signal and background models 
in the (p(top|D), m& t ) plane is performed to extract the top quark mass. 

The event selection is based on that used for the cross section measurement (see 
Sec. |5.3.2| ), with the following differences: 



• Four jets (with > 15 GeV and \r)\ < 2.0) are always required. 

• The aplanarity and Ht cuts are not applied. 

• Events containing high-p^ photons or more than one high-p-p electron are rejected. 

These cuts select 91 events, of which seven are tagged. The dominant background is W+jets, 
and about 20% of the background comes from QCD multijets with a fake lepton. 

The events are analyzed using a kinematic fit to the tt — ► ivbbqq hypothesis. As in the 
CDF lepton+jets mass analysis, the jets are first corrected to match the partons, and the 
fit is then tried for all possible jet configurations consistent with the 6-tagging information. 
For each fit, this gives a fitted mass m^ t and a fit x 2 ■ The jet configuration with the lowest 
X 2 is retained, and events are required to have this minimum x 2 < 10. After this final cut, 
77 events survive, of which five are tagged. 

The kinematic fit is tested on tt Monte Carlo samples both with and without QCD 
evolution and gluon radiation, together with the simulation of detector effects. The mat 
distributions from these tests are shown in Fig. |42|. The mat distribution at the partonic 
level, without QCD evolution, shows a very sharp peak (with ~ 2.4 GeV/c 2 width). About 
80% of the time the configuration with the lowest x 2 is the correct one. The width of 
the peak is mainly due to the widths of W bosons. Like CDF, D0 finds that, after QCD 
evolution and fragmentation effects are included, the width of the met distribution grows to 
~ 26 GeV. (The jets are reconstructed by clustering particles in a cone of width R = 0.5.) 
With these effects included, the correct configuration is selected only about 40% of the time. 
The reduction in the selection efficiency is due to the confusion arising from the splitting 
and merging of jets, the jet finding efficiencies, and increased jet combinatorics. Including 
detector effects does not change the width appreciably (see Fig. f4~2"lc). It can therefore be 
concluded that the primary contribution to the width of the mass distribution is from the 
extra gluon radiation and jet combinatoric effects, rather than from the detector resolution. 

For the cross section analysis, much of the discrimination between signal and background 
is achieved through the use of the variable Ht- However, this variable is not usable for the 
mass measurement because it is highly correlated with the top quark mass. Therefore, 
making significant cuts on Ht reduces the sensitivity to m t . In light of this, D0 embarked 
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Fig. 42. Monte Carlo distributions of fitted mass for herwig it events with nit = 170 GeV/c 2 , 
showing how radiation and detector resolution affect the width of the mass distribution. The 
hatched histograms show the results when the correct jet configuration is used (regardless of 
whether or not it has the lowest x 2 ). The means and widths are from the Gaussian fits shown by 
the dashed lines, (a) Using herwig partons directly, no radiation or detector effects, (b) Using 
final state Monte Carlo particles clustered into cones. Includes the effects of radiation, but not 
detector effects, (c) After full detector simulation and reconstruction. From Ref. B3- 
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on an extensive search for variables that provide good discrimination between signal and 
background and are weakly correlated with mg t . This resulted in the following four variables: 

• x\ = lfi T . 

• X2 = A, the event aplanarity, as defined earlier. 

• £3 = (Ht — E^! tl )/H z , where Ht = J2Et of all selected jets and H z = J2\E Z \ of 
all objects in the event (lepton, neutrino, and the jets), E z being the momentum 
component of the object in the beam direction. This variable measures the centrality 
of the event. 

m x 4 = &Rjj n ■ Ef in /(E T + where ARjf n is the minimum AR of the six pairs of 
four jets and E™ m is the smaller jet Et from the minimum AR pair. This variable 
measures the extent to which the jets are clustered together. 



The distributions of these variables for signal and background are plotted in Fig. |43. 



We see that, on average, the signal events have larger values of the variables than do the 
background events. But, we also note that, while each variable provides some discrimination, 
none is decisive; one might anticipate that direct cuts on the variables may not be the most 
effective way to use them. It proves to be more effective to treat the variables collectively 
and cut on a single multivariate discriminant T>(x) = s(x)/[s(x) + b(x)], where s(x) and b(x) 
are the signal and background densities. Two multivariate methods are used to approximate 
T>(x): (1) a log-likelihood technique, referred to as the low bias (LB) method (due to the 
small correlation of the discriminant with the fitted mass), and (2) a feed-forward neural 
network (NN). In the LB method, D0 parameterizes the ratios Li(xi) = Si(xi)/bi(xi), where 
Si and bi are the signal and background densities for variable x«. Then one computes L = 
Hi and £>lb = L/(l + L) for each event. The weights Wi are adjusted to minimize the 
correlation between the discriminant and the fitted mass. The NN method, by construction, 
takes into account all correlations between the variables used. A three layer feed-forward 
neural network, with four input nodes, five hidden nodes, and one output node is trained 
on samples of Monte Carlo tt (m t = 170 GeV/c 2 ) and background events. The neural 
network is trained using the back-propagation algorithm, choosing the output to be unity 
for the signal and zero for the background. For a given event, the network output T>nn 



directly approximates the ratio s(x)/[s(x) + b(x)]. Figure (£| shows that V LB and £>nn are 
distributed as predicted and provide comparable discrimination. As expected, the signal 
peaks near unity and the background near zero. 

D0 carries out two analyses, based on the 2\b and £>nn discriminants. The data are 
binned in the (Plb/nn, wifit) space. Both analyses use the same binning for m fit (20 bins, 
over the range 80-280 GeV/c 2 ), but differ in how the discriminant is binned. For the £>lb 
analysis, the events are split into two bins, depending on whether or not they pass the "LB 
cut," defined by V LB > 0.43 and H T - E^ 1 > 90 GeV/c 2 . (The latter condition removes 
mostly background; it does not affect the signal much.) For the NN analysis, the events are 
split into ten discriminant bins, as illustrated in Fig. |45|. 

For each method, a top quark mass is extracted from the binned data using a likelihood 
fit to a discrete set of signal and background templates. The likelihood function, based on 
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Fig. 43. D0 distributions of the discriminant variables x%, X2, x%, and X4 (see text for defini- 
tions) for signal (solid histograms) and for background (dashed histograms). All histograms are 
normalized to unity. 
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Fig. 44. Top: D0 distributions of £>lb and 2?nn for signal (hatched, mt = 175 GeV/c 2 ) and 
background (open). The histograms are normalized to unity. Bottom: The distributions of Dlb 
and Pnn for data (histogram), compared to the predicted signal plus background (filled circles) 
and background alone (open triangles). The normalization is from the result of the likelihood fit. 
From Ref. |53|. 
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Fig. 45. Events per bin (oc areas of boxes) in the (^NNj/Wfit) plane for the D0 NN top mass 
analysis, (a) Expected top quark signal (mt = 172 GeV/c 2 ). (b) Expected background, (c) Data. 
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Table 16. Systematic uncertainties in the D0 lepton+jets top quark mass measurement. 



Jet energy scale 4.0 GeV/c 2 
Event modeling 

tt signal 1.9 GeV/c 2 

Background 2.5 GeV/c 2 

Noise/Multiple Interactions 1.3 GeV/c 2 

Monte Carlo statistics 0.85 GeV/c 2 

Fit method 1.3 GeV/c 2 

Total 5.5 GeV/c 2 



Poisson statistics, is derived by a Bayesian method.^ In this fit, the number of background 
events is left unconstrained. Figure ^6| shows the results of the fits and the negative log- 
likelihoods, as functions of the top quark mass. For the LB method, the fit yields m t = 
174.0 ± 5.6 (stat) GeV/c 2 and the number of background events n b = 53.2lg°3 7 (out of 
77). For the NN method, the corresponding results are m t = 171.3 ± 6.0 (stat) GeV/c 2 
and rib = 48.2lg 1 f t . The total systematic uncertainty is estimated to be 5.5 GeV/c 2 (see 
Table 16), of which 4.0 GeV/c 2 comes from the jet energy scale uncertainty and 3.1 GeV/c 2 
from event modeling (including the modeling of initial and final state radiation in tt events). 

D0 combines the LB and NN results, taking into account their mutual correlation of 
88%, to get an overall measurement of m t = 173.3 ± 5.6 (stat) ± 5.5 (syst) GeV/c 2 . For 
D0, this result shows a marked improvement in precision with respect to the measurement 
published at the time of the discovery. This bodes well for the measurement of the top quark 
mass in Run 2 by both collaborations. 



6.4. Summary of Mass Measurements 

CDF observes 8 dilepton events with an estimated background of 1.3 ± 0.3 events, and 
obtains a mass of mt = 167.4 ± 10.3 (stat) ±4.8 (syst) GeV/c 2 . In the lepton+jets channel, 
76 events are observed, of which 15 are singly tagged and 5 doubly tagged with the SVX. The 
mass obtained from this sample is 175.9±4.8 (stat)±4.9 (syst) GeV/c 2 . The 136-event all-jets 
sample, with a background of 108 ± 9 events, yields m t = 186 ± 10 (stat) ± 8 (syst) GeV/c 2 . 
Combining all the CDF measurements yields m t = 175.3 ± 6.4 GeV/c 2 . 

D0 observes 6 dilepton events with a background of 1.4 ± 0.3 events, and measures a 
mass of m t = 168.4 ± 12.3 (stat) ±3.6 (syst) GeV/c 2 . In the lepton±jets channel, 77 events 
are observed (of which 5 are tagged), with an estimated background of 52 ± 9 events. This 
gives m t = 173.3 ±5.6 (stat) ±5.5 (syst) GeV/c 2 . When one combines the D0 dilepton and 
lepton±jets results, one obtains m t = 172.1 ± 5.2 (stat) ±4.9 (syst) GeV/c 2 . 

Combining all five measurements yields^ m t = 173.8 ±5.0 GeV/c 2 . The measurements, 



summarized in Fig. [47], are in striking agreement with each other as well as with the indirect 
measurements from precision electroweak data. It is remarkable that the top quark mass is 
now known to a precision (3%) that is far better than what was anticipated a few years ago! 
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Fig. 47. Summary of all the top quark mass measurements from CDF and D0. 
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Table 17. Data for the CDF \V t b\ analysis from the lepton+jets (W4J) and dilepton (DIL) 



samples. From Ref. 109 





Bin definition 






Sample 




SVX tags 




SLT tags 


W4J 




DIL 


none 




none 


126 




6 


none 




one 


14 






one 






18 




3 


two 






5 









Total 




163 




9 



7. OTHER STUDIES 



7.1. Measurement of \Vtb\ 



t 



The analyses discussed so far have assumed that the branching fraction of the decay 
► Wb is essentially 100%. Clearly, we would like to be able to test this assumption. 
The mixing between the three quark generations is described by the Cabibbo-Kobayashi- 
Maskawa (CKM) matrixS 




V u d V us V u b 
V c d V cs V c b 
Vtd Vts Vtb 




(26) 



where the unprimed letters denote mass eigenstates and the primed ones denote weak eigen- 
states. The ratio R of the branching fraction for a top quark decaying to a 6-quark to that 
of a top quark decaying to any down-type quark q can be written as 



R 



B(t -> Wb) 


v tb 


2 


B(t^> Wq) 


V td 


2 + 


V ts 


2 + 


Vb 


2 



(27) 



(assuming that top quark decays to non-W final states are negligible). In the Standard 
Model with three generations, the CKM matrix Eq. (p6|) must be unitary. With this con- 
straint, current measurements imply that 0.9991 < \V t b\ < 0.9994 (90% confidence level) If , 
however, the assumption of three generations is not made (while maintaining the unitarity 
of the expanded mixing matrix), then the 90% confidence level interval for \Vtb\ opens up 
to 0.05 < \V t b\ < 0.9994. Removing the unitarity assumption leaves \V t b\ unconstrained. 
Therefore, if a measurement of |Vtb| were to yield a result significantly different from unity, 
this would be indicative of new physics beyond the Standard Model (such as the presence 
of a fourth generation). 

The CDF collaboration has made a preliminary measurement of R using their combined 
lepton+jets and dilepton data samplesi^ltOl As shown in Table [17], CDF divides the data 
into several disjoint bins, depending on the observed fe-tags. For each bin i, the number 
of events expected rij can be calculated as a function of the total numbers of signal and 
background events expected in the samples, the 6-tagging efficiencies and fake rates, and R. 
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These expectations are compared to the number of events actually observed in each bin 
Ni using a likelihood function: 

£ = HPiNiln^llGixjlxj^j), (28) 

where P is the Poisson distribution P(N\n) = e~ n n N /N\. The quantities Xj are other 
variables on which the predictions rij depend, each with a measured value Xj and uncer- 
tainty Oj. They are convolved into the likelihood with a Gaussian function G(x\x,a) oc 
exp[(x — x) 2 /2cr 2 ]. A maximum likelihood fit is then performed to extract R and the quanti- 
ties Xj. For this fit, the number of top quarks in the samples is taken to be a free parameter, 
so that the result will not depend on the tt production cross section. 

The resulting likelihood is shown in Fig. |48|. The fit yields R = 0.99 ± 0.29. The 
uncertainty includes both statistical and systematic uncertainties, but is dominated by the 
statistical component. The 90% and 95% confidence limits for R are, respectively, R > 0.64 
and R > 0.58. 

If three-generation unitarity is assumed, then the denominator in Eq. (|27| ) is unity. The 
result then corresponds to |Vt&| = 0.99 ±0.15, or \V t b\ > 0.80 and \V t b\ > 0.76 at the 90% and 
95% confidence levels, respectively. If the unitarity condition is relaxed, no statement can be 
made about \V t b\ without further assumptions. CDF sets \V t d\ = 0.009 and |Vt s | = 0.04, from 
the midpoint of their 90% confidence levels determined with the unitarity assumption. The 
90% and 95% confidence limits on \V t b\ are then \V t b\ > 0.055 and \V t b\ > 0.048, respectively. 

The dominant component of the systematic uncertainty on R is the 6-tagging efficiency. 
This is measured from the data with a precision limited by the available statistics of the 
control samples. Therefore, the total uncertainty on R should decrease as 1/vN as statistics 
are increased. For a 2 fb _1 Run 2 at the Tevatron, the achievable precision for the R 



measurement!!!!! should be about 2%, corresponding to a 95% CL limit of \V t b\ > 0.20. 



An alternate technique which should become possible in Run 2 is to measure \Vtb\ from 



single top pro duct ion. I!!2I This will be discussed further in Sec. |S.3| . We also note that it is 
possible to obtain indirect information on \ V t b\ from the precision measurements which have 
been made in the electroweak sector. A recent reportE^ analyzes the \Vtb\ dependence of the 
single-loop corrections to the Zbb vertex. From a combined analysis using data from LEP, 
SLC, Tevatron, and neutrino scattering experiments, the authors derive \V t b\ = 0.77^o;24, 
independent of unitarity assumptions. 



7.2. Flavor-Changing Neutral Current Decays of the Top Quark 

Flavor-changing neutral current (FCNC) decays can be used to probe for new physics at 
mass scales which are otherwise inaccessible. For example, the absence of the FCNC decay 
K\ — ^ fi + jjT was indicative of the presence of the charm quark, even though the charm quark 



is several times heavier than the kaon.!!! 3 ] For the top quark, the SM rates for FCNC decays 
are very small,^ with the branching fractions B(t — > cZ) and B(t — > cj) ~ 10~ 13 -10~ 12 , 
B(t eg) ~ 10~ 10 , and B(t —>■ cH) ~ 10~ 14 -10~ 13 . Any observation of these decays at 
future experiments would therefore be evidence for new physics beyond the standard model. 
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Because of its large mass, the top quark provides a unique system in which to look for 
FCNC contributions. In many models, FCNC terms are suppressed by ratios of the quark 
masses to a weak interaction mass scale.El Various models have been proposed in which 
the FCNC decays of the to p qua rk are enhanced. These include models with two Higgs 



doublets," supersymmetryj^S and exotic fermions.tll These models typically predict 
increases of the FCNC decay rates of 3-4 orders of magnitude, which are still too small 
to be seen at any experiments planned in the near future. However, in supersymmetric 



models with baryon number nonconservat ion ,1113 B(t — > cry) can be as large as ~ 2 x 10 



which could be observable at the upgraded Tevatron or the Large Hadron Collider, under 
construction at CERN. 

Of particular interest are the decays t — > cry and t — > cZ, as they would yield the most 
distinctive final states.Ela Indirect limits on these branching fractions may be derived from 
lower energy data; the results are Bit —>■ cry) < 1-2 x 1CT 3 and Bit — » cZ) < 0.04. 

The CDF collaboration has searched for these decay modest using their full Run 1 
sample of ~ 110 pb _1 . The branching ratios are measured relative to a normalization sample, 
consisting of the tagged lepton+jets tt candidates. For the t — » qy mode, CDF uses two 
different cut definitions, depending on whether the W boson (from the normally decaying 
top quark) decayed leptonically or hadronically. For the leptonic case, the requirements are 
a central lepton (either an electron or muon) with > 20 GeV/c, $ T > 20 GeV, a photon 
with px > 20 GeV/c, and at least two jets with > 15 GeV. A photon is identified from 
a cluster of energy in the electromagnetic calorimeter with no matching charged track. For 
the hadronic case, CDF requires a photon with px > 50 GeV/c and at least four jets with 
Et > 15 GeV, one of which must contain an SVX 5-tag. In both cases, the invariant mass 
of one jet and the photon is required to lie in the range 140-210 GeV/c 2 , consistent with the 
mass of the top quark. For the hadronic case, the tagged jet must be one of the remaining 
jets, and these remaining jets must satisfy Et > 140 GeV, consistent with the decay of a 
second top quark. Acceptances relative to the SM decay mode are calculated using isajet 
and a parametric simulation of the CDF detector; the leptonic and hadronic cases constitute 
60% and 40% of the acceptance, respectively. The expected backgrounds are less than half 
an event in each channel. 

A single event is observed in the leptonic channel and none in the hadronic channel. 
The event is kinematically consistent with a t — > W~!ry decay, although the photon momen- 
tum (p T = 88 GeV/c) is quite large for this decay. CDF does not subtract the expected 
background, and thus the single observed event gives a branching fraction limit (95% CL) 
of: 



B(t -> wy) + B(t -> cry) < 3.2%. (29) 

The systematic uncertainties in the acceptance calculation are dominated by the statistical 
uncertainty on the number of events in the normalization sample (34 candidates), followed 
by the uncertainty in the 6-tagging efficiency. 

For the t — > qZ mode, CDF searches for the channel in which the Z boson decays 
into a pair of leptons (electrons or muons) and the other top quark decays into three jets. 
The selection requirements are four jets with Et > 20 GeV and \rj\ < 2.4 and a pair 
of opposite-charge, same-flavor leptons with invariant mass in the range 75-105 GeV/c 2 . 
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Acceptances are calculated as before; the total expected background is 1.2 events, with the 
major components being Z + multijet production and ft dilepton decays where the lepton 
pair invariant mass is close to the Z boson mass. 

A single Z — > fifi event is observed, with event kinematics that are more consistent with 
Z + multijets than with an FCNC top quark decay. Again, no background subtraction is 
performed, and the branching fraction limit (95% CL) is 

B(t -> uZ) + B(t -> cZ) < 33%. (30) 

This search is less sensitive than the one for the t —> qj mode because of the small branching 
fraction of Z bosons into charged leptons. The sources of systematic uncertainty are similar. 

These limits can be improved during future Tevatron runs. For a 2 fb _1 Run 2, limits of 
B{t —> qj) < 3 x 10~ 3 and B(t — ► qZ) < 0.02 are probably achievable.tli 



7.3. Top Quark Decays to Charged Higgs 

The minimal Standard Model contains a single complex Higgs doublet, giving rise to a 
single physical neutral Higgs scalar H°. However, the form of the Higgs sector is not presently 
constrained by experiment, so it is interesting to consider models with an expanded Higgs 
sector. In particular, supersymmetric models require two Higgs doublets, one coupling to the 
up-tvpe qu arks and neutrinos, and the other coupling to the down-type quarks and charged 
leptons. cHU In such theories, there are five physical Higgs bosons: two neutral scalars H° and 
h°, a neutral pseudoscalar A , and a pair of charged scalars H . The relevant free parameters 
are the masses of the Higgs bosons and tan /3, the ratio of the vacuum expectation values of 
the neutral components of the two doublets. 

If rriH+ < rn t — m>b, then the decay t — ► H + b can occur, competing with t — > W + b. 



The branching ratio is symmetric in logtan/3 with a minimum at tan/? = ^Jm t /mb. At the 
minimum, B(t — > W + b) ~ 100%, but the situation reverses for tan/3 — > or tan/3 — > oo, 



where B(t — > H + b) rs 100%. See Fig. 49. The principal decays of H + are expected to be 
H + — > cs, dominant at low tan/3, and H + — > t + u, dominant at high tan/3. In addition, a 
recent papeili 2 ! has pointed out that the virtual t quark decay H + — > bbW + is important 



for m,H+ ~ 135 GeV/c 2 and small tan/3. This is also illustrated in Fig. 



It was also recently realized that there are problems with excluding very high and low 
values of tan/3. This is because the coupling of H + to fermions becomes arbitrarily large 
in these limits, causing the perturbative calculations of the branching ratios to break down. 
The widths of the Higgs boson and the top quark also become very large (10-20 GeV/c 2 ) 
in these regions. The tan/3 range over which the coupling is perturbative is approximately 
0.2 < tan/3 < 170. 

Some early limits on the charged Higgs mass are summarized in Ref. |124| . At present, the 
best limit independent of tan/3 is from LEP 2,EiI!m#+ > 54.5 GeV/c 2 (95% CL), extending 
the LEP 1 limits of iriH+ > 44.1 GeV/c 2 . (Preliminary results presented at the 1998 
Moriond conference^ further extend this to > 56.7 GeV/c 2 .) Limits at high tan/3 



have been derived at the SppStHa and from earlier Tevatron runsHH A limit based on an 



independent analysis of CDF data is described in Ref. 130, and an analysis of r lepton 
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Fig. 49. Regions in the (mjy+,tan/3) plane where the branching ratio B(t — ► H + b) > 50%. At 
low tan/3, ff + decays to cs and bbW + , with the latter mode dominating for m H + > 135 GeV/c 2 . 



At high tan/3, .£f + decays to t + v. From Ref. 122. 
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decays which yields a limit on m#+ at large tan f3 is presented in Ref. |131|. In addition, the 



CLEO collaboration, analyzing b — > 37 decays and assuming a two-Higgs doublet extension 
of the Standard Model, derives a limit@ of m H + > [244 + 63/(tan/3) L3 ] GeV/c 2 . However, 
models with additional structure, such as supersymmetry, can evade this limit .EH 

The CDF collaboration performs a direct search^ for t — > H + b in the large tan (3 region 
by looking for H + — > r + z/, considering tt final states where either one or both top quarks 
decay via a Higgs boson. CDF requires one r lepton with p T > 10 GeV/c, two jets with 
E T > 10 GeV, and one or more additional objects, either leptons (including r) or jets, with 
Pt > 10 GeV/c. At least one of the jets must have an SVX 6-tag. If the charged Higgs boson 
mass is close to the top quark mass, however, the 6-jets may fall below the threshold. To 
regain acceptance for this case, CDF also considers a second ("tt") final state, with no jet 
requirements but with two r leptons with pt > 30 GeV/c. In addition, the two r leptons 
must not be opposite in azimuth (A<ft TT < 160°). The r identification algorithm used is 
similar to the calorimeter-based selection described in Sec. [3.3.4| . 

For both final states, it is required that $ T > 30 GeV and A</>/(l°) + $ T /{1 GeV) > 60, 
where A0 is the distance in azimuth between the $ T and the closest other object. This cut 
removes events with mismeasured $ T . In addition, any event containing a e + e~ or fi + n~ 
pair with mass between 75 and 105 GeV/c 2 is removed, in order to suppress background 
from Z boson production. 

In 100 pb -1 of data, CDF finds seven candidates passing the above cuts (six rjjj events, 
one rjje event, and no tt events). The dominant source of background is fake r leptons, 
estimated from jet data to be 5.4 ± 1.5 events for both channels. Additional background 
sources, which can produce real r leptons, include electroweak W/Z + multijet processes 
and diboson production. These backgrounds are estimated using Monte Carlo simulations; 
the results are 1.9 ± 1.3 and 0.08 ± 0.06 events, respectively. The total background for both 
channels is thus 7.4 ± 2.0 events, consistent with observation. 

Acceptances are computed using isajet. The resulting exclusion contour is shown in 
Fig. |5^ (left) for two different assumed top quark production cross sections. The lower of 
the two, on = 5.0 pb, corresponds to the theoretical prediction for a m t = 175 GeV/c 2 
top quark,c3 and the second is 50% larger than that. For large tan/?, CDF excludes a 
charged Higgs boson with m#+ < 147 GeV/c 2 at the 95% confidence level. This limit can 
be extended by requiring that the product a t jB(tt — > W + bW~b) be consistent with the (then 
current) CDF measured cross section0 of 6. 81^4 pb. The result is shown in Fig. ^0] (right). 
Note that these limits may not be reliable for the nonperturbative region tan/? > 170. 

A second analysis strategy is the "indirect" search, which uses the results of the tt cross 
section and top quark mass measurements. It is based on the observation that the standard 
tt selection cuts are less efficient for t — > H + b than they are for t —>■ Wb. Therefore, if 
t — > H + b were to occur at a significant rate, there would be a shortfall in the measured 
cross section relative to the SM calculation. As an example, consider m#+ = 70 GeV/c 2 
and tan/3 = 0.4 (with m t = 175 GeV/c 2 and a t j = 5 pb). With these parameters, one 
has B(t -> H + b) w 100% and also B(H + -> cs) w 100%. Therefore, tt pairs will decay 
nearly always into a six-jet final state, containing neither a high-p^ lepton nor any significant 
]prp. The acceptance of the lepton+jets analyses for these events is practically zero, so all 
observed events must have been background. But the number of events actually observed 
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is significantly above the background expectation; therefore, this point in the (ra#+, tan/3) 
plane can be excluded at a high confidence level. 

Both the CDFil and D0E3 collaborations have carried out such an analysis. The results 



from CDF for the low tan (3 region are shown in Fig. ^Tj. D0 uses this method in both the 
low and high tan /3 regions, as shown in Fig. Note that the D0 results have not yet been 
updated to take into account the three body Higgs decay and the nonperturbative regions. 

CDF performs one other variation on this analysis, which avoids the need to assume 
a it cross section. This is based on the observation that a large B(t — > H + b — ► csb) will 
suppress dilepton events more severely than lepton+jets events. CDF first measures o t - t 
using lepton+jets events, then uses this to find the expected number of dilepton events, 
taking into account the selection efficiency as a function of m#+ and tan/3. The results of 
this analysis are also shown in Fig. ^T]. 



7.4. W Boson Helicity Fraction in Top Quark Decays 

The W bosons from top quark decay can be either transversely or longitudinally polar- 
ized. The Standard Model prediction for the fraction of longitudinally polarized W bosons is 
F = m 2 t /(2M^ + m 2 t ) ta 0.70 ± 0.01 for m t = 173.8 ± 5.0 GeV/c 2 . The charged lepton from 
the decay of a transversely polarized W boson is preferentially antiparallel to the direction of 
the boost from the top quark rest frame to the W boson rest frame. Longitudinally polarized 
W bosons, on the other hand, tend to emit the charged lepton perpendicular to the boost 
direction. The consequence is that the pt spectrum of leptons from transverse W bosons 
will be harder than that from longitudinal W bosons. Thus, one can extract F Q by fitting the 
observed lepton spectrum to the sum of the contributions from longitudinal W bosons, 
transverse W bosons, and background. (The sensitivity of this measurement could conceiv- 
ably be improved by using information from the mass fit. However, that would reduce the 
available statistics and introduce additional biases which would need to be understood.) 
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width is predicted to be larger than 15 GeV/c . From Ref. 123 
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This analysis has been carried out by the CDF collaboration,!!^ which does a simultane- 
ous fit for F in both its e/i and lepton+jets data samples (of seven and 91 events, respec- 
tively). The fit results are shown in Fig. |53|; the extracted F is 0.55±0.32 (stat)±0.12 (syst). 



7.5. Studies of tt Kinematics 

Armed with the results of the kinematic fit from the mass analysis, one can study the 
kinematics of tt events. Typically, one picks the jet configuration which gives the lowest 
X 2 solution and then calculates quantities based on those jet assignments. For exa mp le, 
the tt invariant mass from the CDF collaboration using its tagged lepton+jets sampleli 3 -^ is 
shown in Fig. EM. It should be realized that the quantity plotted is not what one would 



directly calculate for the invariant mass of the two quarks, but is smeared by the effects 



of QCD radiation, jet misassignments, and detector effects. Figure |55] shows a selection of 
kinematic results from the D0 collaboration!! using its lepton+jets samples. In all cases, 
the data are in good agreement with the SM expectations. 

Recently, the CDF collaboration has also presented a measurement of the mass of the 



hadronically-decaying W boson in lepton+jets tt decays.^ Figure pB| shows the distribution 
of the invariant masses of all dijet pairs in lepton+jets events which pass the cut + 
Y.E^ > 310 GeV. The contributions from non-tt background and incorrect combinations 
are subtracted, and the result is then fit to a Gaussian. The result is 77.1 ± 3.8 (stat) ± 
3.6 (syst) GeV/c 2 , and the significance of the excess is 2.8a. They also have 11 events with 
two 6-tagged jets. (When looking for the second tag, the "jet probability" tag algorithm is 
usedJI3 in addition to the usual SVX and SLT tags.) Figure |57| shows the distribution of the 
invariant masses of the pair of untagged jets in each event. When fit to a Gaussian signal 
model plus background, the resulting W boson mass is 78.1 + 4.4 (stat) ±2.9 (syst) GeV/c 2 . 
The significance of the peak is 2.9cr. Finally, these two analyses are combined, giving a 
W boson mass of 77.2 ± 3.5 (stat) ± 2.9 (syst) GeV/c 2 with a total significance of 3.3cr. 



8. PROSPECTS FOR RUN 2 
8.1. Accelerator and Detector Upgrades 

The Tevatron collider Run 1 ended in February, 1996. The accelerator and both detectors 
are presently undergoing major upgrades for the next run, due to start in the year 2000. 

The major improvement to the accelerator complex will be the addition of the Main 
Injector. This 3319 m circumference machine is composed of conventional magnets and has 
a maximum energy of 150 GeV and a minimum cycle time of ~ 1.5 s. It will replace the 
Main Ring, and will serve as the injector into the Tevatron. It will also be used to produce 
antiprotons, thus removing the detector deadtimes incurred by running beam in the Main 
Ring during collider operations. (It can also be used to support fixed target experiments 
running concurrently with the collider.) The Main Injector is expected to increase the peak 
luminosity of the Tevatron from the present value of ~ 2x 10 31 cm~ 2 s _1 to ~ 8x 10 31 cm _2 s _1 . 
A second machine, the Recycler, will be built in the same tunnel as the Main Injector.EH 
It will be constructed almost entirely from permanent magnets, and will operate at a fixed 



102 



35 

30 

25 

20 

15 

10 

5 


14 

12 

10 

8 



Longitudinal W Fraction in Top Decay (F ) 

(CDF Preliminary) 



Lepton + Jet Channel 



+ Data 



— Best Fit 




left handed W 
longitudinal W 



background 



Dilepton Channel 



Combined Result 
F, = 0.55 ±0.32 ±0.12 




20 



40 



60 



80 



100 



120 140 



160 



180 200 



Lepton P, (GeV/c) 
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From Ref. 123 
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top quarks (two entries per event), (c) The difference in pseudorapidity Ar/ between the two top 



quarks, (d) The difference in azimuthal angle A</> between the two top quarks. From Ref. 53 



energy of 8 GeV. Its principal function will be to serve as a high-reliability repository 
for antiprotons. It will also be used to capture unused antiprotons from the Tevatron at 
the end of a store, recool them, and use them again in another collider store. Use of the 
Recycler is expected to increase the peak luminosity to ~ 2 x 10 32 cm~ 2 s _1 . In addition, 
the superconducting magnets of the Tevatron will be operated at a lower temperature; this 
should enable an increase of the center-of-mass energy from yfs = 1.8 TeV to yfs ~ 2.0 TeV. 
This corresponds to about a 40% increase in the it production cross section. The goal for 
the integrated luminosity in the next collider run is 2 fb _1 per experiment. Also envisioned 
is a possible Run 3, with additional improvements to the proton and antiproton sources, 
which could reach an integrated luminosity of ~ 30 fb _1 . 

Both collider detectors will be upgraded^!! for Run 2. The D0 detector will acquire a 
central magnetic field, a silicon vertex detector, and a scintillating fiber tracker, while CDF 
will upgrade to an expanded vertex detector that will provide coverage of the entire luminous 
region of the beam. These improvements will significantly increase the efficiency for tagging 
6-jets in it decays. Many other improvements will also be made to both detectors in outer 
tracking, calorimetry, and muon coverage. 



8.2. Cross Section and Mass Measurements 



CDF and D0 have made impressively precise measurements of the top quark mass: 
m t = 175.3 ± 6.4 GeV/c 2 , and m t = 172.1 ± 7.1 GeV/c 2 , respectively. The corresponding 
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cross sections are a t j = 7.6i};5 pb and 5.6 ± 1.8 pb. These results are based on samples that 
range in size from about five to one hundred events. In the next run, sample sizes should 
increase by at least twenty-fold. 

The po tential for tt physics at the Tevatron during Run 2 and beyond has been studied 
extensivelytlJ^ The expected yields for each detector are on the order of 160 dilepton 
events, 1200 lepton+ > 4 jets events, and 500 double-tagged lepton-l- > 4 jets events. 
For the cross section measurement, the statistical component of the uncertainty scales as 
1 / y/N : . The dominant contributions to the systematic uncertainty are the uncertainties in 
the tt acceptance and background estimates. Both of these are amenable to study using 
control samples from the data, and therefore can be also expected to scale as 1/ \/N . The 
limiting factor could be the error on the luminosity, which can be measured to ~ 5% using 
the W — > Iv rate. With these assumptions, the precision of the cross section measurement 
in Run 2 should be about 8-10%. 

For the mass measurement, the statistical uncertainty should again scale as 1/yN. How- 
ever, the present measurements are already limited by systematic uncertainties, so the chal- 
lenge will be to reduce this component. One of the dominant uncertainties is that in the 
jet energy scale. However, given large statistics, the energy scale can be well characterized 
using Z + multijet events. With an adequate sample of double 6-tagged events, it should 
be possible to use the hadronically-decaying W bosons in tt — > £ + jets channels to perform 
an in situ calibration of the energy scale. Therefore, the uncertainty on the jet scale should 
decrease as 1/yN. The other major uncertainty is that in modeling QCD radiative effects 
in tt decays. Presently, these are modeled using parton-shower Monte Carlo programs such 
as herwig. However, some theoretical progress has recently been made in understand- 
ing the phenomenology of gluon radiation in tt events,0illl and the situation should be 
greatly helped by the availability of tt data with sufficient statistics to provide meaningful 
constraints on models. We therefore expect this uncertainty also to decrease. The total 
achievable uncertainty in the top quark mass measurement from Run 2 should be about 
3-4 GeV/c 2 for each experiment. 

A precise measurement of the top quark mass, along with that of the W boson provides 
a crucial test of the electroweak theory as well as a strong constraint on the mass of the 
Higgs boson. Figure [58] shows the correlation between the W boson mass and top quark 
mass for a range of Higgs boson masses. Also plotted are the current direct and indirect 
measurements of m t and My/. 



8.3. Single Top Quark Production 

So far, this review has focussed on the QCD pair production of top quarks. As mentioned 
in Sec. ||, top quarks can also be produced singly, through the electroweak U W*" (Fig. |J) 
and jy-gluon fusion (Fig. [5]) processes. (There is also a contribution from pp —>■ tW, 
but the rate for that process is very small at the Tevatron.) The cross sections for these 
processes have been calculated at NLO. The results for m t = 175 GeV/c 2 and y/s = 2 TeV 



are a{pp — > tb(tb) +X) — 0. 9 pb for the W* process,!!!^ and a{pp — > tq{tq) + X) = 2.4 pb for 



the VT-gluon fusion process.^ These should be compared to the resummed NLO tt cross 
section of ~ 7 pb.0 
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Fig. 58. The correlation between the W boson mass and the top quark mass as predicted by 
the electroweak theory of the Standard Model, for various possible masses of the yet undiscovered 
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The single top quark production cross section is directly proportional to |V^| 2 . There- 
fore, a measurement of the cross section gives a measurement of |V#,|, independent of any 
assumptions about the unitarity of the CKM matrix. Note that this strict proportionality 
fails in the presence of certain forms of new physics, such as anomalous couplings of the 
top quark. However, evidence for such new physics should also be present in quantities 
which do not depend on \V t b\, such as the ratio of the two subprocess cross sectionsjl'lll! so 
the effects can, in principle, be disentangled. 

One would search for single top quark events by selecting events with exactly one high- 
Pt lepton, large $ Tl and exactly two jets, one of which has a 6-tag. This sort of selection 
would yield a signal/background ratio of about 1 : 2,0 the dominant background being Wbb 
production, with smaller contributions from tt pair production and other QCD processes. 
The signal can be further enhanced by plotting the distribution of m(Wb) and looking for 
an excess at around the known mass of the top quark, as shown in Fig. |59|. 

In Run 2, with 2 fb _1 , 1^1 can probably be measured to a precision of ~ 10%.^ 
(The limiting factor in this measurement is expected to be the error on the total integrated 
luminosity, which has been assumed to be 5%.) 



8.4. Top and New Physics 

The large mass of the top quark is certainly allowed in the Standard Model, but is not 
required. In fact, the Standard Model offers no explanation as to why the top quark is 
so heavy. But because of its large mass, the top quark may be a sensitive probe of new 
physics. The heaviness of the top quark may, in fact, reflect the presence of new physics at 
the electroweak scale. Furthermore, as we discussed in Sec. [3], in supersymmetric theories, 
a heavy top quark can induce the spontaneous breakdown of electroweak symmetry. If this 
were found to be true, it would represent a major step forward in our understanding of the 
origin of mass. We mention below some popular theories of new physics, and outline the 
ways in which its manifestations might be studied using the top quark. 

A very promising class of theories is based on the concept of supersymmetry. (See 



Ref. |146| for a recent review.) There are many valid and intriguing reasons for the prevailing 
optimism concerning their veracity: supersymmetric theories permit the unification of the 
strong, weak, and electromagnetic coupling constants at a grand unification (GUT) scale, 
consistent with the experimental bound on the proton lifetime; they solve the problems of 
gauge hierarchy and naturalness. To add to their allure, superstring theories, which may 
lead to a theory of all forces, predict supersymmetry! 

Supersymmetry (SUSY) is a conjectured symmetry between fermions and bosons that 
predicts the existence of a supersymmetric partner for each particle of the Standard Model. 
(SUSY can also be construed, in a more profound way, as an augmentation of spacetime with 
fermionic dimensions.) In some supersymmetric models, top quarks can be produced in the 
decay of heavy supersymmetric particles; the top quark itself may also decay into lighter su- 
perpartners. The study and precise measurement of top quark production rates can provide 
useful information about such nonstandard processes. It is expected that the production of 
single top quarks, proceeding through electroweak processes, will be particularly sensitive 
to new physics. Supersymmetric theories usually impose conservation of "i?-parity," which 
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Ill 



is a quantum number defined such that it has a value of R = +1 for normal particles and 
R = — 1 for their SUSY partners. However, i?-parity conservation is not required, and if it is 
not imposed, the rate of single top quark production qq — > bt can be greatly enhanced over 
that predicted by the Standard Model. Such reactions would be induced by the exchange 
of a supersymmetric quark (a "squark") in the t-channel, which violates baryon number, 
as well as by the lepton- number violating exchange of a supersymmetric lepton ( "slepton" ) 
in the s-channel. The decay of a top quark into its supersymmetric partner (called a top 
squark or a "stop") and the lightest neutralino, a process predicted in i?-parity violating 
SUSY models, could be observed during the next run of the upgraded Tevatron. 

Another attractive, perhaps radical, approach to understanding electroweak symmetry 
breaking and the origin of fermion masses and mixings, is that of new strong dynamics. 
In these theories, new strong interactions of fermions and bosons are invoked at the scale 
A ~ 1 TeV, and cause the dynamical breaking of electroweak symmetry. Consequently, 
no elementary scalar bosons (such as the Higgs) need exist. The popular models of new 
strong dynamics are those of technicolor, topcolor, and their variants. In the technicolor 
modelJi^J which is analogous to QCD, the chiral symmetry is spontaneously broken through 
a technifermion condensation mechanism, thereby giving masses to the gauge bosons. The 
fermions acquire masses through new interactions (extended technicolor) with the tech- 
nifermion condensate. To account for the heaviness of the top quark, one postulates the 
existence of new topcolor interactions.!^! In these models, many signatures would involve 
the top quark in the final state. Particularly interesting and unique are signatures of res- 
onances that decay to it. In technicolor models]^ a spin-zero color-octet resonance, the 
techni-eta, is produced in gluon-gluon collisions and decays into it or a gluon pair. Such a 
resonance would be seen as a distortion or broadening of the it or two-jet invariant mass 
distributions, but its effect on the bb mass distribution would be negligible. In topcolor 
models,^ a spin-one coloron (top-gluon), produced mainly by quark- ant iquark annihila- 
tions, would decay into it and bb with roughly equal probability, and would appear as a 
broad resonance in both distributions. If an enhancement were to be observed in the it 
mass spectrum, it would be important to study it at the Large Hadron Collider (LHC) 
due for completion in 2005, at CERN. This machine will provide pp collisions with a center- 
of-mass energy of = 14 TeV, with luminosities of up to 10 34 cm" 2 s _1 . At the LHC, 90% 
of the top quark pairs are expected to be produced by gluon-gluon fusion, in contrast with 
the Tevatron, where 90% of the top quark pairs come from quark-antiquark annihilations. 
The aforementioned models have very different predictions of rates and consequences at the 
two colliders. The two colliders hold promise of great discoveries and perhaps some surprises. 

Finally, as Chris Quigg has noted recently,^ in spite of its fleeting existence, the top 
quark may, in fact, have a profound effect on the basic properties of the everyday world. At 
the dawn of the third millennium, with the improved Tevatron and upgraded CDF and D0 
detectors, we expect to refine the measurements of the properties of the top quark and to 
study its dynamics, not only for its own sake, but also in the hope that we shall be afforded 
the good fortune to unravel some of nature's mysteries. 
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